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The structure and thermodynamic state of a system changes under the influence of external electric 
fields. Neutral systems are characterized by their dielectric constant e, while charged ones also by 
their charge distribution. In this Colloquium several phenomena occurring in soft-matter systems 
in spatially uniform and nonuniform fields are surveyed and the role of the conductivity <r and 
the linear or nonlinear dependency of e on composition are identified. Uniform electric fields are 
responsible for elongation of droplets, for destabilization of interfaces between two liquids, and for 
mixing effects in liquid mixtures. Electric fields, when acting on phases with mcsoscopic order, 
also give rise to block copolymer orientation, to destabilization of polymer-polymer interfaces, and 
to order-order phase transitions. The role of linear and nonlinear dependences of e on composition 
will be elucidated in these systems. In addition to the dielectric anisotropy, existence of a finite 
conductivity leads to appearance of large stresses when these systems are subject to external 
fields and usually to a reduction in the voltages required for the instabilities or phase transitions 
to occur. Finally, phase transitions which occur in nonuniform fields are described and emphasis 
on the importance of e and <r is given. 
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I. INTRODUCTION 

There are numerous ways the structure and phase of 
a system can be influenced by external fields. The grav- 
itational field can lead to changes in structure and com- 
position, and in many cases even to phase transitions. 
However, while omnipresent on Earth, its effect is usually 
quite small. A magnetic field has several advantages, but 
it is also weak unless the system is magnetic. Shear forces 
are effective in changing the phase and structure of soft 
materials, but their presence is undesired in many cases. 
Electric fields, on the other hand, have several interest- 
ing properties: they have a strong effect even on neutral 
materials, they can be switched on or off, and they are 
ideally suited for the nanoworld, because an electric field 
increases linearly with a decrease in the electrode size, 
for a given electrode potential. 

The subject of this Colloquium is the influence that 
electric fields have on interfaces between two liquids or 
polymers. This is a broad subject and we do not aim 
at fully covering it but rather mention in this review 
the main effects. An electric field deforms the inter- 
face and may give rise to a dynamical instability. Sim- 
ilarly, fields tend to orient ordered phases of heteroge- 
neous polymeric materials in such a way as to minimize 
the electrostatic energy. Uniform electric fields can also 
lead to the creation of interfaces between liquids. Lastly, 
spatially nonuniform fields have a strong influence on the 
thermodynamic properties of liquid and surfactant mix- 
tures. 

We adopt a continuum coarse-grained approach, where 
all quantities vary smoothly enough to be described by 
continuous fields. In addition, interfaces between two 
pure immiscible liquids are taken to be infinitesimally 
thin. We seek to find the spatial distribution of liquid or 
polymer composition and electric field. The main ingre- 
dients are the dielectric constant difference Ae, e" , and 
the conductivity a, to be defined below. The system's 
behavior depends on the these quantities and also on 
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whether the fields are spatially uniform or not. This ap- 
proach neglects the molecular details and therefore lacks 
accuracy, but it is general enough and captures the phys- 
ical mechanisms at play. 

Let us first recall few basic laws of electrostatics. The 
electric field E is derivable from an electrostatic potential 
ip, such that E = —Vip. The potential satisfies Poisson's 
equation V • (eV"0) = — p, where p is the free charge den- 
sity and e is the medium dielectric constant. We assume 
throughout this article that the system is isotropic, and 
therefore e is not a tensor but a scalar. The displacement 
field D is usually given by the linear relation D = eE and 
is the quantity conjugated to the electric field. 

The component of D perpendicular to a sharp inter- 
face between two materials with dielectric constants £\ 
and £2 is continuous across the interface. Thus, the per- 
pendicular component of E is discontinuous and is given 
by E\E\ — £2-^2, where E\ and E2 are the electric fields 
in the perpendicular direction. In contrast, the parallel 
component of E is continuous across the interface, while 
that of D is not. In systems where the charge on the 
conductors bounding the dielectric material is given, D 
is the natural variable, and E = E(D). The electro- 
static energy density of the pure dielectric is given by 
F cs = (1/2) J (D 2 /£)d 3 r. Conversely, in situations where 
the potential is specified on the set of bounding conduc- 
tors, E is the natural variable, and D = D(E). A Leg- 
endre transform then yields the electrostatic energy, 



£E 2 d 3 r 



(1) 



Note the minus sign in front of the integral. Here is a 
simple argument for it; a full explanation can be found 
in Landau et al. (1984). Our system is subjected to a 



fixed voltage. We can therefore imagine a large condenser 
of capacitance C and charge Q connected to our system 
in parallel. The voltage imposed by this condenser is 
Vo = Q/C. Initially, our small system, having a small 
capacitance c <C C, is not charged (E = everywhere), 
and the total electrostatic energy is U\ — (1/2)Q 2 /C. 
We now connect the two condensers, and charge q <C Q 
enters our system. We find q = Qc/(c + C) and U% = 
(l/2)Q 2 /(c + C) ~ (1/2)Q 2 /C(1 - c/C). To first order 
in c/C, the change in electrostatic energy is found to 
be U2 — U\ = — (1/2)cVq. This argument shows that 
for a system under fixed voltage, the minus sign in Eq. 
([IJ correctly accounts for the work done by the external 
power supply. 

In this article we are interested in liquids, liquid mix- 
tures, and block copolymers in electric fields. It is 
convenient to define an order parameter <p(r), a spa- 
tially dependent dimensionless quantity denoting the rel- 
ative composition of one liquid or copolymer component, 
< 4> < !• We denote by ip the variation in 4> from the 
average value <po, 



0(r) 
Mr)> 



0o + ¥>(r) , 
. 



The variation ip induces a variation in the dielectric con- 
stant £. If ip is small enough, one may write a constitutive 
relation e(0) as a Taylor series expansion to second order 
in <p, 



£(0) = £ + Aep+ ^e'V 



(3) 



£ = e((j>o) and is the average dielectric constant if e" is 
absent from the expansion. At the moment we consider 
"neat" dielectrics, which contain no dissociated ions; 
presence of salt will be allowed later. The electric field 
originates from the presence of a given arbitrary collec- 
tion of conducting bodies at fixed, prescribed, potentials 
and/or charges. We define Eo as the electric field which 
is present in the system when e is constant everywhere, 
£ = £. Variations in composition ip lead to variations in 
£, and since £ and E are coupled via Laplace's equation, 
one has variations in electric field. 

We may thus write to quadratic order in ip 



E — Eo 



1 9 
E x p+ -F, 2 <p 



(4) 



When we later refer to "uniform electric fields" , we mean 
that E is constant everywhere. This means that E 
actually originates from a parallel-plate condenser or 
from nonideal planar electrodes lying far enough from 
the point of interest, such that field inhomogeneities can 
be safely neglected. Of course, even if the zeroth-order 
field Eo is uniform, composition variations lead to field 
nonuniformities, as is evident in the above expansion. 

Let us look at the different terms in an expansion of 
the electrostatic energy density [Eq. dlj] in powers of <p, 



fcs — — -eEq — 



eEq • Ei - 
2A£E • Ei + -e"^ 



AeE 2 1 p 
'E 2 -EqV 



1 

2 

2 , 



(£El 

0(p 3 )- (5) 



(2) 



(Note that energy densities are marked by lowercase let- 
ters.) 

The first term on the right is an unimportant constant. 
The two terms in linear order of ip are inconsequential 
for the thermodynamic state of the system as long as the 
external field Eo is uniform. To see this, one may write 
Ei as a sum of two components: Ei = anEo + aj_Ei_L, 
where Eu_ is the component perpendicular to Eo. If Eo 
is uniform, one finds that a\\ is independent of r and 
J <pEi ■ Eod 3 r = 0. As a results, the linear term in ip in 
Eq. (I5J) vanishes upon spatial integration, recalling that 
(p(r))= 0. When Eo varies in space, this is no longer 
true since a dielectrophoretic force acts on the system. 
Several drastic thermodynamic changes become possible, 
as is discussed in Sec. |V] 

The first and second terms in the second line of Eq. ^ 
(oc p 2 ) are important in uniform fields. The second term 
is twice as large as the first one and opposite in sign, and 
the two sum to give a free energy contribution propor- 
tional to the dielectric contrast squared +(As) 2 . This is 
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a free energy penalty for dielectric interfaces perpendic- 
ular to the external field. We explain how these terms 
give rise to a normal-field instability in liquids (Sec. 
to various orientation effects occurring in mesoscopically 
ordered polymer phases (Sec. |III[ ) , and to the phase be- 
havior of liquid mixtures and block copolymer melts (Sec. 
IV ). The so-called Landau mechanism term, proportional 
to e", is responsible for a modification of the liquid- vapor 
and liquid-liquid coexistences in electric fields (see Sec. 

Tvl 



II. NORMAL FIELD INSTABILITY 

In this section we describe the interfacial instability 
occurring when an initially flat interface separating two 
immiscible liquids or a liquid and a gas is subjected to a 
perpendicular electric field. 



A. Dielectric interfaces 

When an initially spherical liquid droplet is put in a 
uniform external electric field, it elongates in the direc- 
tion of the field. The degree of elongation is given as a 
balance between electrostatic energy, preferring a long, 
needlclikc, drop, and surface tension, preferring a spher- 
ical objec t (|Taylor[[l964l). |Q'Konski and Thacher| ( |T953"1 ) 
and later Allan and Mason ( 1962 ) obtained the following 



expression for small deformations of drops: 
R\\ - R± _ 9 Re 2 E 2 {ex - e 2 ) 2 



Ri 



■R± 16 7 (ei + 2e 2 ) 2 



(6) 



£i and e 2 are the dielectric constants of the drop and 
the embedding medium, respectively. Rn and R± are 
the radii of the ellipsoid in the directions parallel and 
perpendicular to the electric field of strength Eq far from 
the drop, R is the unperturbed drop radius, and 7 is the 
interfacial tension between the two liquids. Clearly, the 
drop elongates in the field's direction. The deformation 
vanishes in the absence of dielectric contrast, that is when 
£1 = e 2 , and does not depend on the sign of E . The case 
of a conducting drop is also given in the limit e±/e 2 — > 00. 

Some of the results of Allan and Mason for dielectric 
drops showed an opposite behavior that could not be ex- 
plained by the theory of neat dielectric liquids - in a few 
cases the drops became oblate rather than prolate. These 
results led G. I. Taylor to propose his "leaky dielectric" 
model. Taylor realized that even a small conductivity 
of the embedding liquid could lead to significant changes 
to the drop shape. Indeed, the tangential stress cause by 
the ionic flow leads to flow inside the drop (Taylor 1966 1. 



The solution of the full clcctrohydrodynamic problem led 
him to suggest a discriminating function $ obeying 



$ = R(D 2 + 1) - 2 + 3(RD - 1) 



2M + 3 
5M + 5 



(7) 



Here, D, R 7 and M are the ratios of dielectric constant, 
resistivity , and v i scosit y of the outer liquid to those of 
the drop (Saville 19971. Drops are prolate when $ > 0, 
oblate when $ < 0, and spherical if $ = 0. Qualita- 
tive agreement has been found between experiments and 
theory derived from Taylor's initial study (Saville 19971. 

The role of residual conductivity, as understood by 
Taylor, will be further highlighted in this review. Let 
us, however, return to pure dielectric fluids. Consider 
this time two incompressible fluids denoted 1 and 2, con- 
fined in a condenser with plate separation L, area S, and 
potential difference V [see Fig. [lja)]. The location of 
the flat interface between the fluids is denoted y = h. 
The permittivities are £1 and e 2 and the electric fields 
are E\y and E 2 y, respectively. Continuity of eE across 
the interface together with hE\ + (L — h)E 2 = V gives 
us Ei = s 2 E /(ei + hAe/L) and E 2 = eiEi/e 2 , where 
Eq = V/L is the average electric field and Ae = e 2 — e\ 
is the "dielectric contrast" . The electrostatic energy per 
volume of the condenser is 



F cs {h) 
LS 



1 £i£ 2 E 2 

2 £1 + hAe/L 



(8) 



As a result, if Ae > 0, the system will reduce its energy 
if h becomes as small as possible, h = 0. The maximum 
energy is when the interface is at h = L. The electrostatic 
pressure on the interface is 



p es (h) 



1 dF cl 



1 £i£ 2 A£ffg 

S~dh~~ ~2{ei + hAe/Lf 



(9) 



The sign reflects the fact that the electrostatic force tends 
to thin the film if Ae > 0. 



(a) 



(b) 



(c) 



■ y=L 

■ y=h 
1 v=0 



FIG. 1 Two liquids in uniform electric field, (a) Schematic il- 
lustration of two liquids with permittivities e\ and £2 confined 
between two plates at separation L and voltage difference V . 
The pressure on the interface located at y = h is given by 
Eq. (b) An initially flat interface at ft = L/2 can break 
into two parts. The electrostatic energy, given by Eq. (10 1, 
prefers h — or h = L. The preferred state is shown in (c). 



Let us now go a step further in elucidating the role 
of dielectric interfaces by allowing the flat interface to 
break into two parts of equal area with heights h and 
L — h, such that the volume of the fluids is conserved 
[Fig. [ljb)]. For simplicity we assume the initial interface 
to be at y — L/2. We ignore edge effects and use the 
previous exercise for the electrostatic energy to get 



F es (h) 
LS 



e\e 2 eE 2 



2e l£2 + (h/L)(l-h/L)(Ae) 



(10) 
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Here, F es is symmetric around h = L/2, and the analy- 
sis shows us that the minimum energy is achieved when 
h = (h = L is equivalent). At this state, the fluid 
interface perpendicular to the field has disappeared and 
was replaced by an interface parallel to the external field 
[Fig. 0Jc)]. The difference between the initial (h = L/2) 
and final (h — 0) states is 



AF CS 
LS 



1 (Ae) 5 



(11) 



This important result shows us that (i) dielectric inter- 
faces parallel to the external field are electrostatically 
favored over perpendicular interfaces and that (ii) the en- 
ergy difference between the two cases scales like Eq and 
is proportional to (Ae) 2 . The deformation of a dielectric 
drop in electric field is in line with this understanding 
- by elongating along the field, the drop decreases the 
area of surface perpendicular to the field compared to 
the unperturbed sphere. 



B. The instability in pure dielectric liquids 



y=0 



>l'=0 




FIG. 2 Stabilization and destabilization by electric field, (a) 
Two liquids in perpendicular electric field, ho is the unper- 
turbed interface location. Interface perturbations are unsta- 
ble if their wavelength is long enough, (b) A liquid film is 
stabilized in tangential electric field. 

The above derivation leads to the intuitive understand- 
ing that there is an electrostatic free energy penalty for 
dielectric interfaces perpendicular to the external field. 
Thus, the perpendicular field destabilizes the fluid inter- 
face initially para l lel to the s ubstrates, as found by|Tay-| 



and Melcher and Smith ( 1969 ) 



lor and McEwan (1965 
[Fig. |2j»J. Melcher a 
trie field stabilizes an interface parallel to it (Melcher 



so noted that an imposed elec- 



and Schwartz 



1968 Onuki 1995b). In the following we 



present a simplified treatment of the dynamical destabi- 
lization process following the lines of |Herminghaus| ( |1999[ ) 



and Schaffer et al. ( 2000 ) . The main aim is to find how 



the period and time constant of the pattern evolving in 
the film depend on the surface tension and electric field. 

Consider fluid confined between two parallel and flat 
electrode separated by a distance L and potential differ- 
ence V, as is schematically depicted in Fig. [2^a). The 
fluid has dielectric constants e\ and unperturbed thick- 
nesses hg. We consider first the thin-film case, where the 



gap contains a gas with dielectric constant e 2 - Gener- 
alization to the bilayer case of two viscous fluids is rela- 
tively easy. As a simple approximation, the flow profile is 
assumed to be Poiseuille-like, and the integrated current 
along the y direction, j(x), satisfies 



j{x) 



h 3 dp 
3ry dx 



(12) 



where rj is the liquid viscosity and x is the direction par- 
allel to the interface. 

There are three different contributions to the pressure 
on the film p: the first one comes from the interfacial 
tension 7 between the two fluids, pi n t = —jh'(x). The 
second contribution is electrostatic, as given by Eq. pj. 
For thin enough films, van der Waals forces come into 
play, giving rise to a disjoining pressure pdis = A/6h 3 , 
where A is the Hamaker constant. For thin liquid films, 
gravity may be neglected; however, gravity effects may 
easily be incorporated (Onuki, 1995b I. 



One may consider perturbations of the flat interface 
h(x) = h + Sh(x, t), where 5h <C fig. The linear stability 
analysis will be restricted to the long wavelength limit, 
Sh'(x) <C 1. We may thus expand p cs and pdis to linear 
order in Sh/ho to obtain 



p(x) ~ —jSh" 



A 

2hJ 



£l e 2 (Ae) 2 E 2 
L (ej + h Ae/Lf 



Sh " 1 " zv — — — - — ^Sh + const. 

(13) 

This expression for the pressure, together with Eq. 



(12 1, is used in the continuity equation dh/dt + dj/dx = 
0. We look at surface waves of the form Sh(x, t) = 
h q e l l r cos(qx), where r is the characteristic exponential 
time for the modulation with wave number q and h q is 
the amplitude. The dispersion relation between r and q 
is readily obtained to be (Schaffer et al. , 2001 1 



3i] 



■ 2 (c 2 



(14) 



The generalized "healing length" is defined by the rela- 
tion 



A 



e 1 e 2 (Ae) 2 E$ 



27^0 jL^ + hnAe/Lf 



(15) 



A positive r means the modulation 5h(x, t) grows in 
time; negative r shows exponential decay. Clearly, all 
q's smaller than are unstable, r is infinite when 
q — > since liquid must then be transported to very long 
distances; the opposite limit, q — > 00, is also reasonable 
because very short wavelengths are rapidly attenuated 
due to surface tension. 

The fastest-growing mode q* is given by q* — ^ x /V2- 
Since A ~ 10~ 20 J, films thicker than few nanometers are 
dominated by the electrostatic forces even at moderate 
field strengths E . If the dispersive part can be neglected, 
we find that the wavelength of the fastest-growing wave 
A* is 



X* 



2n/q* - j 1/2 /AeE 



(16) 
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Thus, the most unstable wavelength can be reduced by 
increasing the dielectric contrast Ae or the field or by 
decreasing the surface tension between the two liquids 7. 

This normal-field instability also occurs in bilayers of 
two viscous liquids. A straightforward generalization 
(Lin et al. 20011 gives the same expression as in Eq. 



(141 



only with the viscosity 7/ replaced by a rather com- 
plicated function of the two liquid viscosities, (7(771,772)- 
The time scale of the dynamics is accelerated by as much 
as 50 times, while q* and A* stay intact. The main ad- 
vantage is the possibility to use liquid pairs with small 
intcrfacial tension, thereby reducing A*. 



C. Experiments 

Initial experiments were conducted on thin polymer 
films above the glass transition temperature, with elec- 
trode gap L < 1 /im. In addition to the basic under- 
standing of interfacial phenomena in electric fields, ex- 
periments were motivated by the possibility of a new 
lithography technique ( |Chou and Zhuang]|1999 1. Indeed, 
amplification of the most unstable mode leads to a film 
with a well-defined periodicity. Figure [3] shows a plot of 




(£,£j/£ 2 



2002) 



FIG. 3 Plot of the most unstable wavelength A* vs reduced 
electric field for various thin films. E\ and E2 are the 
electric fields in the two liquids, Ao = V 2 (Ae) 2 /yje^) 1 ^ 2 , 
and E = V/Xo. The expected relation from Eq. (16 1, 
A*/A = 2tv(E 1 E 2 /E 2 )~ 3/ ' 1 , is the straig ht line with slop e 
of —3/4. Adapted with permission from (Lin et al. 
Copyright 2002 American Chemical society. 



A* for a series of experiments with different parameters: 
liquid thickness ho, spacing L, voltage V, etc. The plot 
shows a good agreement with Eq. (16) when the axes 
are suitably defined (|Lin et al. 2001 1. For a summary of 



data from several groups, the reader is referred to Pease 
and Russell (120031). 



Hierarchical hexagonal structures have been obtained 
in sophisticated experiments by the use of a trilayer, 
namely, a liquid- liquid-air sandwich (Morariu et al 



riodicities, depending on the dynamical process and on 
the relative volumes of the components. 



1. Topographic electrodes 

With smooth electrodes, the pattern period A* [Eq. 
( fl6] )] depends on the electric field, interfacial tension, and 
dielectric contrast. For technologically motivated rea- 
sons, it may be beneficial to achieve different patterns. 
The use of a topographically patterned electrode enables 
facile and rapid duplication of a "master mask" onto the 



liquid film (Schaffer et al. 2000) since the interfacial in- 



stability is amplified in places where the electrode gap is 
smaller. Figure [4] shows an implementation of this idea. 
Subsequent quench freezes the structure. 



uTwinrinn/ 




J J U J U J 




FIG. 4 Pattern transfer. Use of topographically patterned 
electrodes enables transfer of the top electrode onto the poly- 
mer liquid, as is depicted schematically in (a). Tapping- mode 
atomic force microscopy images in (b) and (c) show the result- 
ing topography of the polymer layer. Adapted by permission 
from Macmillan Publishers Ltd: ( jSchaffer et~oT] |2000[ >. 



2003). The resulting structures have two different pe- 
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as 



FIG. 5 Dispersion relation between wave number [scaled 
by (e 2 V 2 /7i 3 ) 1/2 ] and inverse growth time [scaled by 
Vi7L 3 JslV 4 ]. Curves 1-3 are for pure dielectric liquids, 
Eq. (14 1. Curves 4-6 are from the "leaky dielectric" model 



with small dimensionless conductivity S. The conductivity in 
curve 7 is infinite. Reprinted from (Pe ase and Russet) |2002[ ), 
Copyright (2002), with permission from Elsevier. 



2. Residual conductivity 

Work with polymer pairs gives smaller A* because of 
the reduced surface tension, but at the same time the 
dielectric contrast Ae is diminished. How can this defi- 
ciency be overcome? We must return to the role of resid- 
ual conductivity in polymers. Conductivity difference be- 
tween two layered liquids leads to charge accumulations 
in interfacial regions and to additional interfacial stress. 

Pease and Russel (2002) studied the application of Tay- 
lor's leaky dielectric model to this system. Their analy- 



sis for slightly conductive polymer films yielded growth 
exponents and characteristic wave numbers much larger 
than that of perfect dielectrics: even very small conduc- 
tivities decrease the growth exponent r markedly and 
decrease the fastest-growing wavelength A* by a factor of 
2-4. This is shown in Fig. [5] where q and r have been 
appropriately scaled. The dimensionless conductivity S 
is defined by S — ^jaL 3 /e^V 3 , where a is the conduc- 
tivity. The numerical value of S greatly varies and can be 
between 0.1 and 10 16 . Note that according to the model, 
small conductivity has a drastic effect on the dispersion 
relation and that changing from S <C 1 to S > 1 does 
not change much. 



D. Near-critical fluids 

It is in order to briefly mention here near-critical flu- 
ids. The two fluids described above may be a liquid in 
coexistence with its vapor phase or a binary mixture of 
two partially miscible liquids. If one approaches the crit- 
ical point from below, T — > T c , the surface tension and 
the dielectric contrast vanish like 7 ~ 7 (1 — T/T c ) 2y and 
e 2 -£i oc {T c ~T) f:> , where v ~ 0.625 and /? ~ 0.33 arc the 
exponents characterizing the correlation length and the 
liquid-vapor density differe nce (zer o-field magnetization 
in magnetic systems). As Onuki (1995bj) pointed out, 
this implies that A* can be reduced upon approach to T c 



X* ~ 7 1 / 2 /Ae ~ (T c - Tf 



(17) 



The relaxation dynamics of our simple model is influ- 
enced by the vanishing of 7 as T c is approached and is 
different in the short- or long-wavelength limits. One 
should bear in mind, though, that the derivation used 
above breaks down in this limit since we have used the 
assumption of sharp interfaces between the liquids. The 
model is valid as long as the interfacial width between 
the liquids w ~ (T c — T)^ 1 / 2 is much smaller than the 
unstable wavelength: 



E. Other related instabilities 

It is useful to mention severa l related instabiliti e s orig 
mating from electric fields 



Andclman et al. (19861 



considered a Langmuir monolayer of polar molecules. 
These dipoles are constrained to the interface between 
two fluids, typically water and air, pointing in a di- 
rection perpendicular to the interface. Since paral- 
lel dipoles repel each other, their concentration may 
show modulations. On writing the dipole surface den- 
sity as n(x) — n + n q e lqx , they found the electro- 
static energy of the dipoles per unit area to be F cs — 
— (l/2)|g|£ 2 /[£i(£i+£2)]Gpn 2 . Here d is the electric dipole 
moment of individual molecules. The electrostatic en- 
ergy thus prefers a modulated state with infinitesimal 
period. The system is unconstrained in the direction per- 
pendicular to the two liquids, hence the — \q\ dependence. 
The competition between the long-range dipole-dipole re- 
pulsion and interfacial tension between "phases" , scal- 
ing as q 2 , leads to an unstable mode with finite wave- 
length ( |McConnell et al.\ [1984 1. The simple analyti- 
cal model, valid in the long- wavelength limit, enabled 
them to construct a phase diagram of modulated phases 



with hexagonal and stripe symmetrie s (|Andelman et al. 
1987|. Similar behavior was found by Garel and Doniach 



( 1982 ) for thin uniaxial ferromagnetic slabs subject to 



a perpendicular magnetic field even though the physi- 
cal origin of the dipoles is different in the two cases. 
Another system of interest consists of a charged end- 
group polymer brush placed inside a parallel-plate con- 
denser. The electrostatic energy per unit area associated 
with end-group height undulation of the form h(x) — 
ho + h q cos(qx) can be written to quadratic order in h q 
as F es = — <r 2 /e\q\ cosh(qh ) cosh[q(L — h )]/ smh(qL)h q . 
Here a is the charge per unit area, L is the surface sep- 
aration and £ is the dielectric constant of the uniform 
embedding medium. One retrieves the F QS ~ — \q\ de- 
pendence in the symmetric case where h = L/2 and in 
the limit L — ► 00 ( Tsori et al. . 2008 1 . Again, the surface 



is unstable with respect to a finite wavelength. This re- 
sult essentially recapitulates the studies on the instability 
appearing in 3 He- 4 He interface when it is charged with 
ions ptzet all 119841 iWanner and Leidererl 11979) . The 
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instability of Sec. [TT]is dynamical because once the inter- 
face is parallel to the field, the field stabilizes it. Here, 
however, the instability appears in equilibrium because 
due to its charge, the interface is frustrated even if it is 
perpendicular to the substrate. 



In a related work, Du and Srolovitz (20041 consid- 



ered the stability of the metal electrodes themselves to 
surface modulations. Two parallel metals with poten- 
tial difference V and separated by an insulator are com- 
monly found in many situations, such as in micro-electro- 
mechanical-systems, microswitches, and close to scan- 
ning transmission microscopy tips. They used the same 
long-wavelength approximation and obtained the electro- 
static energy F cs = —(l/8)\q\eEQ coth(qL)h^, where h q 
is now the amplitude of electrode surface modulation. 
When interfacial tension was added, they found that all 
surface modes with q < q c are unstable, where q c is given 
by q c — zEq coth(q c L) /4-f. However, the system is most 
unstable with respect to long wavelengths, that is, q* = 0. 

The stability of a poorly conducting lipid membrane in 
aqueous environments in perpendicular electric field was 
studied as well ( Sens and Isambert . 2002 ) . The accumu- 
lation of charges on the opposite sides of the membrane 
is unbalanced if the membrane is not completely flat. 
The destabilizing electric field acts like a negative sur- 
face tension, tending to enlarge the membrane area. For 
a freely suspended membrane, a hydrodynamic theory 
gives the fastest-growing q mode q* ~ (E R) 2 / 3 , where 
Eq is the external field and R ~ 10 5 is the ratio between 
the membrane and solvent resistivities. The correspond- 
ing wavelength is A* = 2n/q* ~ 0.5/im and the growth 
rate is r ~ 10 5 - 10 6 s" 1 . 

Lastly, we mention the normal-field instability in fer- 
rofluids, an important phe nomenon discovered by |Cow- 
ley and Rosensweig (19671. When the interface between 



a ferrofluid and a nonmagnetic fluid is subject to a per- 
pendicular magnetic field, the surface becomes unstable 
if the magnetization exceeds the critical value, M c . The 
static pattern is hexa gonal, and its period A* is given by 



A* = 2n('j/gAp) 1 / 2 (Andelman and Rosensweig 
Rosensweig 1997[ ). 
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F. Immiscible liquids in electric fields: Electrowetting 

Until now, we have described (i) how an initially flat 
liquid layer becomes unstable under the influence of a 
perpendicular electric field and (ii) the deformation of a 
liquid drop in external fields. We now briefly describe 
an "intermediate" situation, that of a liquid drop placed 
on a solid flat substrate in electric field. Classical elec- 
trowetting describes the change in the wetting properties 
of two immiscible liquids due to the field (Mugele and 



Baret 


2005 


applications 


et al.\ 


2004) 



This is a broad topic with many important 



( Stone 



and steady-state situation. The elongation of dielectric 
drops not in contact with any substrate, as considered by 
O'Konski and Thacher [see Eq. (pi], may lead us to think 
that drop will elongate in the field's direction, thereby re- 
ducing the contact area with the substrate and increasing 
the contact angle. The general shape change is correct, 
but at the contact line this intuition fails: the contact 
angle 9 stays the same. The reason that 9 is independent 
of E Q is because the electrostatic energy scales as the vol- 
ume while the interfacial energies scale as the area. Upon 
looking at ever smaller regions close to the three-phase 
contact line, we thus find that the electrostatic force be- 
comes negligible compared to the interfacial forces. The 
apparent contact angle, measured at a macroscopic scale, 
may be larger than the field-free angle. 

For a conducting drop, the situation is very different: 
charge accumulation at a very thin layer at the substrate 
means the energy contribution of the electric field be- 
comes proportional to the surface and not to the volume. 
This leads to a modification of the interfacial properties 
and hence the apparent wetting angle changes from 9 to 
9*. The main differences from Taylor's work on weakly 
conducting drops are (i) the existence of a solid surface 
and the localization of electric field at it and (ii) the em- 
bedding medium, usually a vapor, is nonconducting. 

Lippmann's original work on electrolytes asserted that 
the solid-vapor interfacial tension 7 SV is unaffected by 
the potential, but solid-liquid interfa cial tension 7 s i i s 
reduced by a value proportional to V 2 (Lippmann 1875). 



This reduction is due to the spontaneous creation of an 
electric double layer at the substrate. However, in the 
cases where the electrode is metal, electrolysis usually 
makes it difficult to achieve high voltages. It is therefore 
beneficial to cover the electrode with a dielectric insulator 
of thickness d and permittivity e. In this situation, the 
dielectric is responsible for the system's capacitance, and 
the effective solid-liquid interfacial tension 7*] is given by 
Berge| ( [19931 ), 

7 S * - 7s. - £V ■ (18) 

The apparent contact angle 9* is found from substitu- 
tion of 7*j in Young's equation and is given by 

e(V~V ) 2 
2jd 



cos 9* = cos 9 - 



(19) 



Consider first a dielectric drop embedded in a di- 
electric medium. For simplicity, we assume dc voltage 



The Vo term accounts for spontaneous charging - it is 
common that a solid in contact with an electrolyte in- 
herits a net charge by ion adsorption or ionization of co- 
valently bound groups. For instance, common glass near 
water ionizes to make SiO - and releases a proton. The 
above supposes a wedge-shaped interface and is correct 
on a macroscopic scale. On a scale smaller than thickness 
of the dielectric layer d, the curvature of the liquid- vapor 
interface is not fixed. 

In the absence of a disjoining pressure, the interface 
shape is governed by the generalized Laplace's equation, 



7«(r) -pes(r) = const. 



(20) 
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Here k is the local curvature and p ca cx V 2 is the electro- 



static pressure, calculated globally for the drop (Mugele 
At the mesoscopic scale, a numer- 



and Baret 



2005) 



ical procedure assuming a circular contact line showed 
that the electric field diverges weakly at the three-phase 
contact line, and the slope approaches Young's angle, 6, 



at the contact line (Bienia et al. 2006 Buehrle et al. 



2003). This theoretical finding has been verified experi- 
mentally recently (Mugele and Buehrle 20071. Note that 



the drop's shape needs not stay circular, and a static in- 
stability of the contact line was observed in high voltages 
( |Mugele and Herminghaus 2002). In contrast to the in- 
stability discussed throughout Sec. [TTJ, here the complex 
interfacial shapes are static due to a balance between 
Laplace and electrostatic pressures. 

The influence of time- varying ac fields depends on the 
field's frequency to: for a liquid drop of dielectric constant 
e and conductivity a, the drop behaves as in dc field if 
i t is in the quasistatic regime, that is, when lu <C lu c 
(Landau et al. 1984 Mugele and Baret 20051, where 



= cr/e . 



(21) 



In the opposite regime, lu 3> u> c , the mobile ions do 
not have enough time to response to the field before it 
changes sign. Electric double layer is not created, and 
therefore the field acts throughout the whole drop vol- 
ume. The electric field thus exerts a body volume force 
again just like for static dc fields in pure dielectrics. Ex- 
perimentally, demineralized water have a ~ 4 x 10~ 6 S 
m _1 and therefore at frequencies ui > 10 4 s _1 the be- 
havior is similar to the dielectric case. The frequency lu c 
appears also as an important measure of the influence of 
ions on orientation of block copolymers by electric fields 
(see Sec. |IIle| . 

As a final comment, we stress that all of the above 
is not valid close to a critical point: if T — T c is small 
enough, the correlation length diverges, and the inter- 
facial tensions 7, 7 SV , and 7 s i become dependent on the 
field distribution and droplet shape. 



III. BLOCK-COPOLYMER ORIENTATION IN UNIFORM 
ELECTRIC FIELDS 

The preceding section described several examples of 
deformation of an interface between two immiscible liq- 
uids or between a liquid in coexistence with its vapor. In 
both cases the two phases are macroscopic, separated by 
a thin interfacial layer. An interesting question regards 
the effect of electric field on complex materials which self- 
assemble into ordered structures with typical lengths on 
the mesoscopic scale. We will concentrate on one such 
material - block copolymers (BCPs). These consist of 
two or more chemically distinct polymer species con- 
nected by a covalent bond. Their self-assembly results 
from a competition between enthalpic interactions and 
chain stretching of entropic origin (IBates and Fredrick- 



research in the past years because many of their interest- 
ing mechanical, electrical, rheological, and other proper- 
ties can easily be fine tuned for optimum performance in 
nanotcchnological applications ( Park et al. 2003 1 . 



The phase behavior of diblock copolymers, made up of 
two polymers A and B, is governed by two parameters: 
</>o, the volume fraction of the A polymer (0 < 4>q < 1) 
and N\- X is the so-called Flory parameter and is ap- 
proximately inversely proportional to the temperature 
and TV is the polymerization index. In symmetric di- 
block copolymers, the A and B parts of the polymer 
have the same length, and 4>o — 1/2. If the temper- 
ature is high enough (small Nx), the melt is found in 
a disordered phase. Cooling below the so-called order- 
disorder temperature (ODT) leads to the appearance of 
a lamellar phase. This is a one-dimensional and periodic 
phase of A- and B-rich domains of width d comparable 
to the polymer radius of gyration. At fixed temperature, 
further increase in the length asymmetry between the A 
and B blocks (larger |0 O — 1/2|) increases the spontaneous 
curvature: the polymers adopt a configuration where the 
shorter block is confined to the inside of cylinders ar- 
ranged on a hexagonal lattice. Even larger increase in 
|0o — 1/2| leads to the formation of spheres arranged on 
a bec phase. Other more complicated ordered phases 
appear in block copolymers and even in simple diblocks 
(e. g., gyroid phase), but we ignore them here. The re- 
gion in the phase diagram in the (cj> ,Nx) plane close 
to the ODT point, given by Nx — 10.5, is called the 
weak-segregation regime. In this region, domain spac- 
ing d scales like TV 1 / 2 . The region farther below the 
ODT point, Nx S> 10.5, is called the strong-segregation 
regime, where the A and B blocks are highly segregated 
from each other, and d ~ TV 2 / 3 . 



A. Orientation mechanism 

The orientation of block copolymers and indeed many 
other soft-matter phases can be influenced by several ex- 
ternal fields. Examples are shear flow, van der Waals or 
chemical interaction with surfaces, confinement by walls, 
etc. At first inspection, it seems that electric fields are 
very weak since the electrostatic energy stored in a molec- 
ular volume vq, v q sE 2 , is much smaller than the thermal 
energy fc^T. However, experiments have shown that this 
estimate is too simplistic and vq is not the relevant vol- 
ume. The first experiments on block copolymers in elec- 
tric fields were carried out by Amundson et al. (19911. 



In a series of papers they developed the key concepts 
for alignment and stability of lamellar phases in elec- 
tric fields (Amundson et al. 1994 19931. In the first 



son 1999 ) . Block copolymers have attracted considerable 



set of experiments, polystyrene-polymcthylmcthacrylate 
(PS/PMMA) was anionically synthesized. Samples were 
heated above the ODT point, held for about 10 min, 
and then cooled below the ODT in the weak-segregation 
region. In the absence of an electric field, a small-angle 
x-ray scattering (SAXS) reveals a characteristic ring, cor- 
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responding to a repeat period of d = 23 nm. Samples un- 
der electric field, however, exhibit strong anisotropy, as 
is manifest in the clear peaks appearing in the direction 
perpendicular to the electric field. 
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FIG. 6 Orientation of PS/P MMA block copolymer lam ellae, 
(a) The bulk experiments of Amundson et al. ( 1993 1 show 



near-perfect alignment of lamellae in the direction parallel to 
the external field, (b) Time evolution of the orientational or- 
der parameter S for two samples with different cooling rates 
as measured by birefringence experiments. Solid line is an 
exponential fit. Adapted with permission from (Amundson 
et al.\ 19931. Copyright 2002 American Chemical society, (c) 



TEM closeup of cylindrical diblock copolymers near an elec- 
trode. Inset is schematic illustration of the full setup. Dashed 
arrows are field lines and box indicates area of image. Note 
how cylinders curve and follow the field lines near the corner. 



From (Morkved et al. 
A A AS. 



1996 1. Reprinted with permission from 



Figure [6][a) is a transmission electron micrograph 
(TEM) of lamellae oriented in electric field. In order 
to quantify the alignment of ordered phases, |Amund-] 



son et al. ( 1993| ) used the quantity S defined as fol- 
lows. Denote by <j)(r) the local A-polymer composition 
(0 < 4> < !)• I n the disordered phase, <f) is equal to its 
average value 4>q. One can write </>(r) as a sum of plane 
waves of the form 



0o + <£>(r) , 

q#0 



«qr 



(22) 



For example, close to the ODT point, a lamellar phase 
can be described by 4>(r) — O + <£l cos(qx), where ipz, 
is an amplitude and the normal to the lamellae is chosen 
in the x direction. Similarly, hexagonal and cubic phases 
contain three and six modes, respectively. For a lamellar 
phase, the average orientation S is then defined by 



In nematic liquid crystals, S is commonly used to quan- 
tify order. In the beginning of the process, there are 
many random grains in different directions not correlated 
with the field, the sample is macroscopically disordered, 
and therefore S — 0. As orientation of lamellar grains 
proceeds, q turns until it becomes perpendicular to Eo 
for reasons which will become clear below. Perfect ori- 
entation has S = —1/2. Figure |6jb) is a time-evolution 
plot of S(t). The evolution is slow because the sam- 
ples are very viscous. The sample with slow cooling rate 
equilibrates faster to the electrostatically preferred ori- 
entation. Figure |6jc) is a TEM image taken by Morkved 
et al. (19961 close to the electrodes. Clearly, lamellae 



seem to follow the curved field lines. 



B. Experiments 

The recent decade has witnessed an explosion of work 
on electric-field effects in block copolymers. Experiments 
can be roughly divided to two types: those carried out 
in a block copolymer melt and those in a solution. In 
the first case, pure block copolymer melt is subject to 
an external electric field. Since the melt is very vis- 
cous, the sample is heated to elevated temperatures and 
subsequently annealed under electric field. SAXS, TEM, 
and small-angle neutron scattering arc used to character- 
ize the copolymer structures. Annealing times are very 
long, and therefore experiments in copolymer melts in- 
trinsically probe static phases. On the other hand, when 
copolymers are dissolved in a solvent, the viscosity is 
not high and orientation kinetics in electric field can be 
recorded. 



1. Statics 



The experiments of Amundson et al. (1 19941) showed 



that block copolymer structures can be oriented in the 
bulk. This work has stimulated further attempts to ori- 
ent thin film structures. The advantage in using thin 
films is their possible technological applications and the 
reduction in required voltage following the reduction in 
size. The disadvantage in size reduction is that the elec- 
tric field, acting throughout the sample volume, needs 
to overcome increasingly more dominant interfacial en- 
ergies acting between the copolymers and the bounding 
surfaces. 

The experiment in Fig.[6jc) and others (Mansky et al. 



(23) 



19981 showed that BCP domains can be oriented later- 
ally along an external field even in a thin film. Consider 
a BCP melt sandwiched between two flat and parallel 
surfaces at distance L from each other and potential dif- 
ference V. One would naively expect that there is one 
transition field, below which BCP lamellae (or cylinders) 
are parallel to the substrate and above which domains 
are perpendicular to it. The electric field acts through- 
out the sample volume, and therefore the order of mag- 
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nitude of field strength required for effective orientation 
is given by eE 2 ~ 'j/L, where 7 and e are the typical dif- 
ference between the interfacial tensions of the polymers 
with the substrate and dielectric constant. In most ex- 
perimental systems 7 ~ 10-100mN, L ~ 10-100 /im, and 
e ~ 5e , and we therefore find E ~ 10 3 -10 7 V/m. The 
larger value of these estimates is quite a large field, but 
it is nonetheless considerably smaller than the dielectric 
breakdown threshold, E ~ 10 8 V/m. 
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FIG. 7 Summary of sample orientation with respect to the 
field. Orientation angle is the angle between cylinders' 
axis and external field for three different film thicknesses L. 
Shaded area is a coexistence "mixed" region, where both ori- 
entations are observed. Critical fields E\ and E2 mark the 
transition fields between the three possible states (cylinders 
parallel or perpendicular to the surface and a mixed state). 



Adapted with permission from ( Thurn-Albrecht et al. 2000 1 



Copyright 2002 American Chemical society. 



vent. The use of a polymer solvent allows additional 
experimental flexibility: for example, a selective solvent 
for one of the polymers may, in principle, effectively in- 
crease the dielectric contrast Ae, and it also enables use 
of high molecular weight and branched polymers. In ad- 
dition, because the solution is less viscous than the melt, 
the full orientation kinetics can be recorded in real time. 
However, a usual deficiency of this approach is the non- 
selectivity of common solvents: as the solvent dissolves 
in the respective polymer blocks, the effective dielectric 
contrast Ae between domains is diminished. 

Bulk samples of polystyrene-poly(2hydroxyethyl 
methacrylate)-polymethylmethacrylate 
(PS/PHEMA/PMMA) with number-average molec- 
ular weight M„ of 82 000 g/mol were dissolved in 
chloroform and put in a cylindrical capacitor with 
average electric field Eq = 1.8 V//im. The solvent 
evaporated in a controlled manner and under the influ- 
ence of external field, and the resulting samples were 
characterized by SAXS, TEM, and differential scanning 
calorimetry. These techniques show that the PHEMA 
block is miscible with the PMMA block and that the 
triblock copolymer actually behaves very similar to 
a PS/PMMA diblock of enhanced dielectric contrast. 
The experiments showed that block-copolymer solutions 
can be efficiently oriented parallel to the external field, 
overcoming interfacial interactions. In addition, they 
also opened a window for quantitative consideration of 
defect statistics and, more importantly, the dynamics of 
grain orientation and defect annihilation. 



Subsequent work on BCP thin films verified that 
electric field can indeed orient lamellae and cylinders 
from a direction parallel to the substrate to a direction 



perpendicular to it (and parallel to the field) (Thurn- 



Albrecht et al. 2000 1 . Anionically synthesized samples 
of PS/PMMA in the disordered (high temperature) state 
were exposed to electric fields of 1-25 V//im. One alu- 
minum electrode was in direct contact with the copoly- 
mers, while the second one was insulated from the copoly- 
mers by a Kapton sheet. The sample was then cooled 
down to the cylindrical phase and held for 14 h, and fi- 
nally cooled to room temperature. SAXS experiments, 
shown on Fig. [7] revealed the following general picture: 
surprisingly, there are two critical fields E\ and E 2 . If 
E < Ei, cylinders are parallel to the substrate, while 
above E\ they are destabilized. If E\ < E < E2, there 
is a mixed state, showing characteristics of both parallel 
and perpendicular cylinders. If E > E2, the sample is 
fully oriented perpendicular to the substrate and parallel 
to the external field. 

As mentioned above, due to the high polymer viscos- 
ity, experiments in the melt can be performed with poly- 
mers of limited molecular weight. Typically, large electric 
fields are required in order to overcome defects that per- 
sist in the sample and to accelerate the slow dynamical 
process. A different approach was presented by Bokcr 



2. Dynamics 

As is explained above, one of the advantages of orienta- 
tion in polymer solutions is the possibility to track orien- 
tation dynamics. In order to record dynamics, a high-flux 
scattering source is required because the exposure time is 
limited as compared to static experiments. Another im- 
portant consideration is the copolymer volume fraction 
in solution: less copolymer not only reduces the viscosity 
but also reduces the effective dielectric contrast Ae. The 
electric field needs to be high enough so that orientation 
terminates before solvent evaporation completes and the 
structure becomes virtually immobile. B5ker et a?.| (2003 1 
used lamellar-forming polystyrene-polyisoprene (PS/PI) 
block copolymer dissolved in toluene. The samples were 
exposed to synchrotron SAXS of high beam energy and 
photon flux during exposure to electric field. The syn- 
chrotron beam direction was perpendicular to the electric 
field's direction. The measured average sample orienta- 
tion, as indicated by S(t) of Eq. (23), showed a similar 
qualitative behavior as in Fig. [6]Jb^ But there is a big 
difference - the dynamics are 20-fold faster. The relax- 
ation of S(t) is characterized by a single exponential with 
time constant r. As expected, r increases with increasing 
copolymer content and decreasing temperature. 



et al. (20021, whereby copolymers are dissolved in a sol- 



Using in situ SAXS experiments, Boker et al. (2003) 
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identified two mechanisms for orientation of ordered 
phases: close to the ODT point, the polymer solution un- 
dergoes a transition between only two orientations (from 
perpendicular to parallel to Eo), with no intermediate 
orientations. Grains of lamellae in the favorable direc- 
tion grow on the expense of unfavorable grains, resulting 
in grain boundary migration. However, far from the ODT 
point (lower temperatures) , the scattering pattern has all 
the intermediate grain orientations and thus reflects con- 
tinuous rotation of grains to the favorable direction. The 
optimum grain orientation is not fully achieved since the 
driving force for orientation, torque, is negligibly small 
at long times. 

What is the nature of the "driving force"? Figure [8] 
shows the collection of results from numerous systems 
with different copolymer solutions. For all four copoly- 
mer solutions, the inverse exponential relaxation time 
t -1 scales linearly with -Eg' 7 , with scaling prefactor de- 
pending on the specific solution. The exponent, being 2.7 
and not 2, is attributed to alignment process governed 
by an energy barrier; this is a feasible explanation in the 
weak-segregation regime, where grain orientation is dom- 

~ 2007) . 



inated by nucleation and growth (Schmidt et al 



A simple scaling form exists validating the importance of 
the dielectric contrast Ae: when the data are presented 
for t]/t against (Ae) 2 ^/^, where r/ is the solution vis- 
cosity, all data except one collapse onto a universal curve. 




^^£ 2 [kV 2 /mm 2 ] 



FIG. 8 Inverse exponential relaxation time r\jr against Eq 
from four copolymer solutions as obtained by SAXS. When 
prop erly scaled, most data sets fall on a single universal curve. 
FromlSchmidt et al. ( 2007 1 . Reproduced by permission of The 



Royal Society of Chemistry. 

A general theory of orientation of ordered phases is 
presented below based on a coarse-grained simplified ap- 
proach. 



C. Theory 

The underlying physical mechanism of copolymer ori- 
entation, sometimes referred to as the "dielectric in- 
terfaces," is described below. This mechanism relies 



on dielectric contrast between the different material do- 
mains. The electric field favors one sample orientation 
over another, the energy difference being proportional to 
(Ae) 2 Eq. Intuitively, we may say that there is a free 
energy penalty for having dielectric interfaces perpendic- 
ular to the field's direction. This is true in both the weak- 
and strong-segregation regimes, but the analytical theory 
is different in the two cases. 



1. Electrostatics of strongly segregated lamellae 



— 1(a) 



I" 





FIG. 9 Capacitor analogy for three competing lamellar struc- 
tures confined in a condenser, (a) Parallel and (b) perpen- 
dicular layering and (c) mixed morphology, consisting of two 
parallel layers of thickness a at the electrodes and perpendic- 
ular lamellae throughout the rest of the film. The equivalent 
electrical circuits are illustrated on the left. 

All interesting effects reviewed here are due to dielec- 
tric anisotropy. One obvious anisotropy in e is due to the 
mesoscopic order, in which e is a function of r. Another 
contribution may come from the inherent anisotropy of 
e in the presence of electric field. Indeed, since polymers 
are large objects with many chemically differing regions, 
e in general is a tensor, and therefore the dielectric re- 
sponse is anisotropic in nature. Here we ignore the ten- 
sorial nature of e since it is not observed in most experi- 
ments. The interested reader may turn to the theoretical 
description of block copolymers with anisotropic e given 
by IGurovichl (fl994l (19951) . 



Strongly segregated lamellae are characterized by a 
square-wave composition profile cj). For the electrostatic 
calculation it is therefore reasonable to take lamellae as 
being composed of pure A and pure B polymers, with 
dielectric constants £\ and £2, respectively, separated by 
a sharp interface. The capacitor model allows the cal- 
culation of electrostatic energy of competing orientations 



(Pereira and Williams 1999 Tsori and Andclman 2002). 



Consider a stack of lamellae parallel to the electrodes 
separated by distance L and potential difference V, as 
illustrated in Fig. [9ja) . Such parallel lamellae may be 
stabilized by preferential interactions with the substrate. 
For symmetric lamellae, we use the continuity of the dis- 
placement field D = £E across dielectric interfaces and 
the total potential drop and find the electrostatic energy 
per unit area to be Fes = ~\C\\V 2 , where the capaci- 
tance per unit area Cm is given by 



Cii 



£l£ 2 1 



(24) 
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Here e = {e\ + £2)/2 is the average dielectric constant. 
Similarly, for the perpendicular arrangement [Fig. K3[b)] 
' ~ V 2 and 



(a) 



we have = -\C 



(25) 



Since C± is larger than G\\, the electrostatic energy 
-F cs is always lower for the perpendicular lamellae, with 
the difference between the two orientations being pro- 
portional to (Ae) 2 . Which of the two phases, parallel 
or perpendicular lamellae, is thermodynamically stable 
that is a matter of a competition between the interfa- 
cial energies, scaling like the substrate area, against the 
electrostatic energy, scaling like the sample volume. It is 
clear that there is a critical voltage above which perpen- 
dicular lamellae will be preferred over parallel ones. This 
voltage depends on Ae, L, and the interfacial energies. 

The two structures above are not the only conceivable 
ones. More complex phases, such as the one depicted 
in Fig. |9jc) may be possible. This mixed morphology 
presents an interesting compromise: the system keeps 
few parallel layers at the substrates, thus optimizing in- 
terfacial interactions, but has lamellae oriented in the 
field's direction in the rest of space, thus minimizing elec- 
trostatic energy throughout most of the film. However, 
there is an energetic penalty per unit area of the film, 
7 T , for the creation of a "T-junction" defect. 

The capacitor model allows us to calculate F cs for the 
mixed state as well. The equivalent electrical circuit is 
that of two "parallel" capacitors and one "perpendicular" 
one connected in series. Ignoring edge or "fringe" effects, 
we may write the capacitance of the mixed morphology 

by 



1 

Cm 



2a 1 



L-2a 1 
L Cl 



(26) 
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FIG. 10 Stability diagrams for confined lamellar phases in the 
strong segregation regime in electric fields, (a) L = lOd and 
is fixed. Small value of Eq leads to parallel layers. Increase in 
Eo leads to a direct transition to perpendicular layers when 
^ — (7bs ~ 7as)/7t i s sma U enough. If 8 is large, 8 > 5*, 
increase in the field induces first a transition to a mixed phase 
and subsequently to a perpendicular orientation, (b) Diagram 
in the Eq-L plane with fixed 8 = 5. Eg is scaled by (7 T /d) 1|/2 . 



Adapted with permission from (Tsori and Andelman 20021. 



Copyright 2002 American Chemical society. 



part (b), the transition field E\, above which parallel 
lamellae become unstable, has undulations. These are 
due to the incommensurability between the lamellar pe- 
riod and the film thickness, causing polymer stretching or 
compression. E\ decays like L™ 1 / 2 - this is easily under- 
stood since the interfacial energy is proportional to the 
area, while the electrostatic energy scales as area xLEq. 
E 2 is the transition field between mixed and perpendic- 
ular lamellae and is virtually independent of film thick- 
ness L. Note that |7 BS — j AS | may be so large that E 2 
is above the threshold for dielectric breakdown (~ 150 
V/fjxa). In these circumstances perfect perpendicular 
lamellae throughout the whole film are not possible. 



The factors in front of the inverse capacitances on the 
right hand side are necessary since the parallel layers 
only occupy a region of width a (equal to several lamellar 
widths). It is easy to show that Cm < C m < C±. 

The free energies per unit area of the three competing 
phases are given by 





F 1 - 



Lf P + 2 7a 



-CnL 2 E 2 



Lf p + 2 lAS +2 lT - l -C m L 2 El 



(27) 
(28) 
(29) 



f p is the polymer free energy per unit volume and in- 
cludes bending and compression. 7 AS and 7 BS are the A 
and B interfacial energies with the substrate. 

The minimal model above allows the calculation of 
phase diagrams as a function of Eq, j bs — j AS , and L. 
Figure [10] shows two such cuts in the three-dimensional 
phase space. The transition fields E\, E 2 and £3 are 
lines where the free energies are equal to each other. In 



Experiments validate the existence of a mixed mor- 
phology. Fig. [TTJshows lamellae sandwiched in a thin film 
in the presence of a perpendicular electric field. Clearly, 
the film's center has lamellae oriented parallel to the 
field's direction, while close to the interface with the sub- 
strate or air, the copolymer component with lowest in- 
terfacial energy is preferentially adsorbed, leading to the 
creation of few layers parallel to the substrate. The defect 
separating the two regions, while not perfectly similar to 
the ideal one depicted in Fig. [9j can nonetheless be as- 
sociated with an energy r y T per unit area so the model 
should still be valid. 

Muthukumar has considered other confined ordered 
phases in electric field by using a similar model. While 
the dependence of the transition field between parallel 
and perpendicular orientations on L is weak for lamel- 
lae, they find that for cylinders the L dependence of 
the energy is more complicated, and the transition field 
changes markedly as the film thickness changes (Ashok 
etoOl200T|. 
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2004b) 



FIG. 11 TEM cross section of thin-film PS/PMMA block 
copolymer annealed in electric field Eg ~ 40 V/^rn. Few 
lamellae lie parallel to the substrate (bottom) and polymer-air 
interface (top), while the rest of the film ha s lamellae paralle l 
to the field. Adapted with permission from ( Xu et 
Copyright 2002 American Chemical society. 



2. Electrostatics of weakly segregated structures 

The above derivations are valid in materials (not nec- 
essarily block copolymers) where the dielectric interfaces 
are sharp and Ae is large compared to the average per- 
mittivity e. Close to the critical point, however, the 
composition differences between coexisting phases be- 
come small. The expressions for the electric field and 
electrostatic energy given below only assume small e(r) 
variations and therefore are general and not restricted to 
BCPs. The smallness of e variations, as compared to e, 
enables us to exactly solve Laplace's equation for an arbi- 
trary spatial distribution e(r). Close to the critical point, 
we may write the order parameter as a sum of Fourier 
harmonics [Eq. (22 1], with ip(r) <C 1. Consequently, only 
linear terms are retained in the constitutive relation e(</>) 
and in the expansion of E in powers of ip [Eqs. ^ and 
Q]. The deviation of the electric field from its average, 
Eitp in Eq. Q, is found by writing Laplace's equation 
V • (eE) = to linear order in <p and is given by ( Amund- 
son et a?.[[l993| , 



(30) 



E r 



Aetpq 



(q • E ) q 



The electrostatic energy per unit volume can be approx- 
imated by 



(Asf 
2e 



q 



(q ■ E )V q 9s_ q + const. (31) 



The constant on the right-hand side is the electrostatic 
energy for a uniform phase with <^(r) = <pQ. Any plane- 
wave composition deviation c/3 q e iq ' r from the isotropic 
background leads to an electrostatic penalty. This en- 
ergy, quadratic in Eq as usual, is also quadratic in y? q . 
This is because in uniform electric fields, composition 
fluctuations ip give rise to dielectric deviations Se cx ip 



and electric-field deviations SE cx ip, so the leading order 
term in the free energy density — (l/2)eE 2 scales as ip 2 . 

The torque N acting on a sample of volume V in ex- 
ternal field Eo consistent with the energy formula [Eq. 
plj)] is fTsori et q?.[ |2003b| : 



N diclcc = -2V^|^T( 



,E ) 



(32) 



r (0> E o) = 5^ q y- q (q-E )qx e 



An important feature of Eqs. ( |31[ ) and (32) is the depen- 
dence on the angle between q and E . If q and E Q are 
perpendicular to each other (dielectric interfaces parallel 
to Eo), the energy is minimal and the torque vanishes. 
When q and Eo are parallel (interfaces perpendicular to 
Eq), the torque vanishes, but the energy is highest. 



Equation (31 1 is a useful formula for the calculation 



of electrostatic effects in soft-matter systems close to a 
critical point. When it complements a proper Ginzburg- 
Landau expansion of the energy in powers of tp, it en- 
ables rather easy derivation of various thermodynamical 
expressions on the mean-field level (Onuki 20021. This 
issue will be discussed in Sec. IIVI 



3. Numerical calculations 

A numerical approach for the calculation of electro- 
static effects in spatially nonuniform materials has the 
advantage that it does not assume that Ae/e is small 
or large. Numerical calculations therefore bridge the 
gap between weak- and strong-segregation theories out- 
lined above. Electrostatic calculations for block copoly- 
mers have been so far implemented by mean-field self- 
consistent field theory (SCFT). In SCFT of block copoly- 
mers, the Hamiltonian contains chain stretching penalty 
and enthalpic interactions, originating from unfavorable 
contacts between chemically different monomers. The 
partition function is expressed in terms of the density op- 
erators of the different types of monomers. The so-called 
restricted chain partition function obeys a modified dif- 
fusion equation (Matsen and Schick 1994). In order to 



account for electrostatic effects, a constitutive relation 
e(<f>) is chosen (usually a linear relation), and the elec- 
trostatic energy F cs in Eq. ([I]) is added to the Hamil- 
tonian. The full set of equations must be solved self- 
consistently together with Laplace's equation and sub- 
jected to the proper boundary conditions on the elec- 
trodes, i.e., Dirichlet, Neumann, or mixed boundary con- 
ditions (|Lin et al.\ 120051 iMatsenl |2006b|). The main ad- 



vantage of the formalism for this study is the exact calcu- 
lation of electrostatic energy - the numerical procedure 
does not rely on the analytical approximations such as 
Eq. pi |. 



While static calculation pertains to equilibrium mor- 
phologies, a certain variant of SCFT allows description 
of the dynamical relaxation toward equilibrium. In the 
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scheme developed (Sevink et al. 19991, the free energy 
is calculated similarly to the static case. The time vari- 
ation in the order para meter <j> is the n assumed to obey 
the following dynamics (Onuki, 2002 1: 



dt 



(33) 



where the chemical potential fi is a functional derivative 
of the total free energy including F cs , M is a proper On- 
sager mobility coefficient (usually taken to be constant), 
and 77 is a noise term satisfying the fluctuation-dissipation 
theorem. This approach does not fully take into account 
the polymer viscoelastic flow and reptation effects. In 
addition, the Amundson-Helfand quadratic electrostatic 



energy approxi mation [Eq. ( 31 1] was used instead of the 



full expression (Schmidt et al. 2007). Still, a reason 



able match with experiment was found: the simulations 
report an exponential relaxation of the order parameter 



S(t) [ Eq. (23 1], consistent with experiments (Boker et al. 
20031. 



It should be stressed that the approach utilized above, 
that e is a simple scalar quantity, may fail with some 
polymers. Recent experiments show a peculiar phenom- 
ena, whereby the lamellar period changes by as much as 
5% unde r the influence of an electric field parallel to the 
lamellae (ISchmidt et al. 2008 I. The researchers attribute 



the thinning of lamellae to an anisotropic response of the 
polymers (PS/PI). In essence, the single chain conforma- 
tions of the PI block are affected by the field, and the 
coil becomes more prolate in a direction parallel to the 



field dGurovich[|1994 l. 

Lastly, the cylinders in Fig. [6jc) follow the curved elec- 
tric field lines. The electric field introduces electrostatic 
(Maxwell) stress in the anisotropic sample, and this stress 
must be balanced by elastic forces. The competition be- 
tween electrostatic and elastic forces has not received 
enough attention in the literature and should be further 
investigated. 



D. Instability of block-copolymer interface in perpendicular 
electric field 

In light of the normal-field instabilities described in 
Sec. [n] one may ask whether a similar interfacial insta- 
bility can also occur in block copolymers. The difference 
between self-assembled materials, such as block copoly- 
mers, and simple liquids is that besides their high viscos- 
ity they feature elastic behavior. Indeed lamellar phases, 
for example, have both bending and compression moduli. 
As a result, any deformation of the interface between the 
copolymers results in deviation from the optimum bal- 
ance between interfacial and entropic energies. It turns 
out that an interfacial instability does occur in BCPs in 
electric field, and the nature of the instability depends 
on the distance from the critical point. 

(1995) studied the strong- 



Onuki and Fukuda 



segregation regime by writing the copolymer composi- 
tion as 4>(x,z) — (f>o + ipLCos[(2w/d)x + u(y)], with tp^ 



as an amplitude and d as the lamellar period, and ex- 
panding the energy in small u. They found an interfacial 
instability with two characteristics: (i) two adjacent A/B 
polymer interfaces are parallel to each other and (ii) the 
most unstable wavelength is long, corresponding to the 
lateral film size [Fig. 12 [a)] . This instability is similar to 



the Helfrich-Hurault instability common in smcctic and 



cholesteric liquid crystals ( de Gennes and Prost 1993 1 



(a) 



(b) 




FIG. 12 Schematic illustration of two normal-field instabil- 
ities in lamellar block copolymers (not to scale), d is the 
unperturbed lamellar thickness and A* is the undulation pe- 
riod, (a) In the strong-segregation regime, adjacent lamellae 
bend in phase with each other. A* is macroscopic, A* ^> d. 
(b) Close to the ODT point, lamellae are easier to deform, 
and adjacent lamellae undulate out of phase with each other. 
A* is equal to the lamellar thickness d. 

A different scenario occurs close to the ODT (critical) 
point. In this regime, the composition variations are very 
weak, and lamellae are easy to deform. A stability anal- 
ysis of the Ginzburg-Landau-like free energy expansion, 
appropriate close to a critical point and on the mean- 
field level, reveals that the most unstable wavelength A* 



is the same as the bulk period d [see Fig.[l2[b)] (Tsori and 



Andelman 2002 1 . This is very different from the insta- 
bility in the liquid case. Moreover, adjacent lamellae are 
out of phase with each other. At the onset of instability, 
the BCP morphology is essentially a superposition of two 
lamellar phases with normals parallel and perpendicular 
to the substrate. These modulations were also corrob- 



orated by a more accurate SCFT study (Matsen 2005 



2006c I , which describes undulatory "peristaltic" modes 
similar to that of Fig. [12] These theoretical predictions 



were essentially verified in experiment (Xu et al. 2004a I 



E. Role of residual ions 

The preceding discussion of block copolymers assumes 
that the polymers are purely dielectric materials and ig- 
nores the possible existence of dissociated ions. How- 
ever, almost all polymers have small conductivity due 
to the finite amount of ions in them. For example, in 
anionic polymerization the reaction typically starts with 
an organomctallic reagent, such as butyl lithium. At the 
end of the reaction, each chain is neutralized with wa- 
ter, and one metal hydroxide ion (e.g., LiOH) is released. 
Thus, the average ion density is one ion per polymer 
chain, which is equivalent to a number density of 10 25 
m~ 3 (~ 0.016M) or charge density of 1.6 x 10 6 C m -3 . 
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Even without electric fields, a salt content adds to the 
molecular polarizability of the polymer, and therefore 
changes the interaction parameter x between polymer 
molecules. In addition, like any chemical impurity, salt 
can change the interfacial interactions of the polymers 
with the substrate. The effect of ions on the BCP phase 



diagram was studied quantitatively ( Epps et al. 2002 
2008[). For a given phase, we are inter- 



2003 



Wang et al. 



ested in the effect ions have on the orientation process in 
external fields. Most of the salt content crystallizes and 
is found in neutral clusters inside the polymer, and only 
a small fraction (10~ 3 — 10~ 6 ) is dissociated and partici- 
pates in the orientation mechanism. Of those ions, some 
are complexed to the polymer to which they have a pref- 
erential solubility, usually the more polar block. They 
increase the effective dielectric constant of the polymer 
and therefore increase the effective Ae. 

The rest of the ions are mobile; when subjected to elec- 
tric field, they give rise to screening. Two extreme cases 
follow: in the first, the screening length A is much larger 
than the lamellar period d and in fact may even be as 
large as the film thickness. When this happens, the ion 
cloud is uniformly distributed throughout the film in the 
more polar polymer, and the system behavior is similar 
to that of neat BCPs, only with effective e's. The critical 
field for transition from parallel to perpendicular lamel- 
lae still scales as E ~ L -1 / 2 . Orientation experiments 
were carried out on PS /PMMA samples doped with ions 
flWang et al] |2006| |Xu et al] |2004a[). Ions reside pref- 



erentially in the PMMA block, resulting in an increased 
Ae and correspondingly smaller fields for the orientation 
process. How much the dielectric constant of the PMMA 
block is enhanced depends on the frequency of applied 

[2009 1. 



field and on the temperature (Kohn et al 



Typ- 
ically, the dielectric constant increases with decreasing 
frequency. Large doping with ions can lead to even a 
fivefold increase in the static dielectric constant. 

The second extreme case corresponds to screening 
length A smaller than d and occurs when salt is added. In 
this case, most of the ions concentrate at the electrodes. 
Here one expects the most polar polymer to wet the sur- 
face due to the preferential solubility of ions. The effect 
of electric field is thus reminiscent of electrowetting and 
can be approximately described by an effective interfacial 
energy between the polar block and the substrate. Note 
that in contrast to the case with neat BCPs, increase in 
the electric field increases the tendency of lamellae wet- 
ting parallel to the substrate. The effective surface ten- 
sion depends on the total ion density and vanishes with 
vanishing electrode potential V. The mixed structure 
depicted in Fig. [9] is possible for intermediate screening 
lengths A > d and depends on the vicinity to the ODT 
point. 

The physics of the problem is rich when the external 
field varies sinusoidally in time. When a time- varying 
electric field is applied, the positive charge will tend to 
migrate parallel to the field, while negative charge will 
move in the opposite direction in an oscillatory motion. 



3 , one obtains the 
11 C/m V s. Thus, 



There are four length scales in the problem: the BCP 
domain size is d ~ 50-100 nm, comparable to the polymer 
radius of gyration. The second is the electromagnetic 
wavelength, 2ttc/ '(u)y/s), where uj is the field's frequency. 
The third length is the distance a single ion drifts in 
one half period of the electric field: TrefiEo/uj, where ft 
is the average mobility. Lithium mobility in PMMA is 
estimated to be ft ~ (4-5) x 10 5 m 2 /J s. Typically the 
wavelength of light is much larger than both the drift 
length and the BCP period. Lastly, the system size L is 
usually the largest length. There are three energy scales: 
one is the energy stored in the dielectric material per unit 
volume, t/dicicc = SEq. The second energy is the thermal 
energy fc^T. The third energy is the amount of heat 
dissipated per unit volume due to Joule heating in one 
field cycle: Uj ou \ e = 2ttEq<t/uj. Using an average density 
of dissociated Li ions n ~ 10 21 m 
average conductivity a = e 2 nft ~ 10 
Joule heating be gins to be important at low frequencies, 
uj < uj c [Eq. (21 1], where here ui c ~ 5-50 s _1 . 

The copolymers are anisotropic materials, and there- 
fore the mobility of ions is spatially nonuniform. Ions 
will consequently drift in the direction which is "easy" 
for them, and this direction may not be collinear with 
the external field. As a result, a torque is applied on the 
sample. The total torque works to orient the sample in 
such a direction that it vanishes and the mechanical en- 
ergy is at minimum. Clearly, the torque depends on the 
external field's frequency: small uj and long period means 
ions have a long way to drift before the field changes sign, 
and hence the torque exerted on the polymer matrix is 
large. At large values of uj, ions are practically immobile 
and contribute nothing to orientation. A detailed calcu- 
lation shows that the torque exerted by moving ions is 
dominant over the "dielectric interfaces" if a /uj 3> e . An 
expressions for the torque in this regime exists, and it 



reads as (Tsori et al. 2003b) 



i 'mns — 



8a v fAfi 



r(0,E o ) 



(34) 



This formula should be contrasted with Eq. (32 1. It is 



valid close to the critical point and under the assumptions 
that positive and negative charges have equal mobilities 
and using a linear dependence of mobility on composi- 
tion: fi((f>) = ft + ipAfi. In this regime, occurring when 
uj < uj c , the mobile ions become increasingly important 
and the concept of dielectric contrast should be replaced 
by mobility contrast. In this regime the sample is being 
heated; in most experimental systems, however, this en- 
ergy dissipation does not present a problem since heat is 
quickly removed from the system. 

The stresses introduced by moving ions and by di- 
electric interfaces (Maxwell stress) can lead to real 
symmetry-changing phase transitions and not just to ori- 
entation. This will be discussed in subsequent sections. 
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IV. CRITICAL EFFECTS IN POLYMER AND LIQUID 
MIXTURES IN UNIFORM ELECTRIC FIELDS 

Up to this point, we described several instabilities and 
orientation effects occurring when two immiscible liquids 
or ordered phases of block copolymers are subjected to 
an external electric field. One may ask even a more basic 
question: how does an electric field affect the relative 
thermodynamic stability of the possible system phases? 
We turn to investigate this question on the mean-field 
level as is the practice throughout this text. 



A. Landau theory and experiments 

We consider a bistable system, able of having two dif- 
ferent states at a certain range of the parameters. For 
concreteness one may think of a binary mixture of two 
liquids. The treatment of liquid- vapor coexistence of 
pure molecules is quite similar, the main difference be- 
ing the finite gas compressibility which means that the 
Clausius-Mossotti relation connects between e and den- 



sity ( Hegseth and Amara 2004 Zimmerli et al. 1999a|b I. 
What is the effect of a uniform electric field on the phase 
di agram of liquid mixtures? T his issue was treated first 
by Landau and Lifshitz ( 1957 ) . For simplicity of exposi- 
tion, we consider a symmetric liquid mixture. The free 
energy density in the absence of field, f m , can be written 
as a Landau expansion valid close to the critical point, 



k B T c fm 2 



IT — T r 



T r 



4! 



(35) 



Here ip is the deviation of the composition from the crit- 
ical composition (1/2), vo is a volume derived from the 
second virial coefficient, and u is positive. We use the 
constitutive relation [Eq. pf] to write the electrostatic 
energy density / es ~ -\E c E$-\£L£yE%-\e"ip 2 Z%. The 
constant term is unimportant and the linear term can be 
eliminated by a redefinition of the chemical potential. 
The only important term is the quadratic one. Inspec- 
tion of / = f m + / cs shows that since / es is small, close 
enough to T c , that / os simply serves to redefine the crit- 
ical temperature. Thus 



T r 



AT r = 



T c + AT C 
v a s"E 2 



2k f 



(36) 



T c and the whole binodal curve are increased if e" is pos- 
itive, and in this case electric field leads to demixing. If 
e" < 0, then T c is decreased and mixing is favored. We 



use the value vq 



10 



-28 



m and average electric field 



Eq ~ 10 V//Ltm to obtain the typical shift to T c due 
uniform fields in either case to be v^e" E^/kg ~ 30 mK. 
Subsequent theoretical work employed a renormalization- 



group (Onuki 1995a I and other approaches (Sengers 



et al. 1980 ) to study the vicinity of the critical point, but 



the main result remained: T c changes by few milikelvins 
under the influence of moderate-to-strong fields. 



The theory should be compared with experiments 
which followed. The first experiment was by |Debye and 



Kleboth (1965) on mixture of isooctane and nitroben- 
zene, whose dielectric contrast is Ae = 32.3eo- Under an 
electric field of 4.5 V//im, they found that T c is reduced 
by 15 mK. These results were later verified with great ac- 
curacy by Orzechowski (199 9). D ebye's experiments were 
followed by a work of Bcagle hole] ( |1981 [) . He worked on a 
cyclohexane/ aniline mixture with Ae = 5.8 and in a field 
of 0.3 V /jitm. T c was found to be reduced by as much as 
80 mK. Early ( 1992[ ) worked on the same mixture with 
similar fields. He found that T c does not change at all and 
hinted that previous re sults were due to spurious heating 
occurring in dc fields. Wirtz and Fuller ( 1993 1 worked 



on n-hexane/nitroethane mixture with Ae = 17.7eo and 
found a reduction in T c by 20 mK. 

The experiments above are in contradiction to theory 
since mixing is observed instead of demixing even though 
e" > 0. In addition, there were wrong signs in some ex- 
periments as well as theory. More importantly, in both 
theory and experiment the change to T c is much smaller 
than 100 mK. The only exception is the work of Reich and 
), who worked on a polymer-polymer mix- 
ture having a lower critical solution temperature (LCST) 
and found that T c changes by about 4 K or even more 
under an electric field of ~ 10 V//xm. Two comments on 
the big effect observed are the following: (i) The large 
molecular volume of polymers means their entropy is re- 
duced. v should be replaced by Nvq in Eq. (36). (ii) 



Heating, if present in a LCST system in the homogeneous 
phase, may lead to demixing, as is observed in these ex- 
periments. 

A possible explanation for the mixing observed in ex- 
periments is due to the dielectric anisotropy not ac- 
counted for by Eq. (36 1. Indeed, just as with block 



copolymers or any other anisotropic material, composi- 
tion variations ip lead to dielectric constant variations 8e. 
Compliance to Laplace's equations means that the elec- 
tric field has variations too, and in linear order SE on —Ss. 



As is clear from Eqs. (22 1 and (31 ), the electrostatic free 
energy contribution / cs scales quadratically with ip. We 
therefore conclude that the prefactor of p 2 in a Landau 
series expansion of f m + / es around the critical compo- 
sition has two terms: one proportional to — e" and the 
second to +(Ae) 2 . A positive e" favors demixing while 
(Ae) 2 term promotes mixing and is usually dominant. 

Light-induced phase transitions in mixtures 

We briefly mention several cases where light influences 
the morphology of mixtures. Light in the UV range was 
applied on polymer blends of stilbene-labeled polystyrene 
and polyvinylmethylether (PSS/PVME). This caused 
phase separation of the blend because trans— >cis tran- 
sition changes both the volume and the enthalpic in- 



teraction between the polymers (Nishioka et al. 20001. 
In addition, photochemical reaction was induced in 
anthracene-labeled polystyrene-PVME blends irradiated 
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by periodic UV light. The distribution of length scales 
in the resulting spinodal decomposition depends sensi- 



of </? q , given in Fourier space by 



tively on the applied frequency (Tran-Cong-Miyata et al. 
20041). 



Even in the absence of chemical activity, an intense 
laser beam can lead to a mixing or demixing phase tran- 
sition. Poly(iV-isopropylacrylamide) (PNIPAM) - wa- 
ter and PNIPAM - D2O mixtures were illuminated by 
tightly focused laser light in the infrared wavelength 



(Hofkens et all 1997 Ishikawa et al. 19961. These mix- 



tures have LCST (inverted) phase diagrams and the ini 
tial temperature corresponded to a homogeneous state. 
Aggregation of the polymer was observed in the region 
illuminated by the laser. Two possible explanations 
were suggested. First, water adsorption in the infrared 
regime causes heating, and in a LCST system this may 
cause phase separation. Second, radiation pressure sup- 
posedly works in conjunction with the heating to sepa- 
rate the polymer from the solvent and to accumulate it 
in micron-scale aggregates. Subsequent work (Borowicz 
et al. 1997) was similar, but now poly(-/V-vinylcarbazole) 



was dissolved in two different organic solvents, cyclohex- 
anone and A^A^-dimethylformamide. This work showed 
that the segregation of polymer from the solvent occurs in 
the absence of heating and is due to the strong laser elec- 
tric field. The underlying physical mechanisms responsi- 
ble for the observat ions were theoretically examined by 
Delville et aT1 ( |1999[ ). In Sec. |V]we give a possible mech- 
anism which may be relevant to equilibrium situations. 



B. Block-copolymer phase transitions 



In Sec. |III[ block copolymer orientation in electric 
fields was described as a result of a "dielectric interfaces" , 
stemming from a term proportional to (Ae) 2 in the free 
energy. Equation (32 1 shows that a BCP sample expe- 



riences a force as long as there are dielectric interfaces 
perpendicular to the external field. Phases with uniaxial 
symmetry thus turn until their axis of symmetry is par- 
allel to Eq. When this orientation terminates, their total 
electrostatic energy is exactly equal to — ^eEq x V, where 
V is the total volume. But there are ordered phases, such 
as the bcc phase of spheres, where frustration always oc- 
curs. As E is increased, the crystal deforms, and the 
spheres elongate along the external field. This is remi- 
niscent of the liquid droplet elongation described in Sec. 



II. A but in addition to the interfacial tension and elec- 
trostatic energies, in BCPs there is also elastic energy 
that needs to be taken into account. With increasing 
value of Eq, the crystal deforms, until, at a certain criti- 
cal value of Eq, its energy is no longer the lowest among 



all phases, and a phase transition occurs (Tsori et al. 
2003a|. 



The phase diagram close to the critical point can be 
calculated by a Landau expansion of the energy in series 



knT c 



IE 



T-T Ci v (Ae)' 



eknTr 



(q-E ) 2 



C(q 2 - ql) 2 } + 0(<4) 



(37) 



The last term in the square brackets is specific to mod- 
ulated phases such as block copolymers, and is absent 
in simple liquids. It expresses the penalty of having q 
modes with q different than the preferred inverse period- 
icity, qo = lit I d. The prefactor C contains information 
on the polymer, such as the radius of gyration. The elec- 
tric field (second term) introduces a penalty proportional 
to the cosine squared of the angle between field and q vec- 
tor. As is mentioned above, this field-dependent term is 
usually quite small: when we replace s" by (Ae) 2 /e in 
the estimate after Eq. ( 36 1 , we get the same millikelvin- 
sized effect of renormalization of T c . For block copoly- 
mers, however, the molecular weight can be large and 
thus, replacing vq by Nvq, one may get a larger change 
to T c , similar to the experiments on polymer mixtures 



(Reich and Gordon 19791 




(b) 



Eo 1 


Dis JV Hexagonal 


E, r 


// \ x n 
' / \ N ^ 


0.5 


' / \ 
II \ v 




/ I \ 

/ distorted \ s \ 




j bcc \ \ 


n 





0.38 0.39 0.4 0.98 1 1.02 1.04 1.06 

<|> K/Xtr 

FIG. 13 Block-copolymer phase diagrams, (a) Phase diagram 
for diblock copolymers in the (4>o,x) plane at fixed electric 
field. Two outer solid lines are zero-field hexagonal-to-bcc 
and bcc-to-disorder phase boundaries. Lines A and B are 
transitions from distorted bcc to the hexagonal and disor- 
dered phases, respectively, and line C is the boundary be- 
tween hexagonal and disordered phases in electric field given 
by (£ Nvo/k B T) 1/2 E = 0.2. They m eet at the triple point 
(0o,A r x)=(O-39,11.43). Adapted from ( |Lin et al|p005| . (b) 
Phase diagram in the {Eo,x) plane and fixed composition 
4>o — 0.3. For low values of field, a distorted bcc phase 
lies between the disordered ("Dis") and hexagonal phases. 
The phase lines terminate at a triple point (xtr, -Etr). When 
So > E t i, the distorted bcc does not exist. Solid lines are from 
analytical one-mode approximation, whereas dashed lines are 
from self-consistent field theory. The axes of each calculation 
are scaled by the relevant values of xtr and Eti . Adapted with 
permission from ( Lin et aE] 20051 and (Tsori et al. 20061. 
Copyright 2002 American Chemical society. 



Figure 13 'a) is a phase diagram in the (</>o,x) plane 
(X ~ l/Tjfor fixed electric field, obtained with the SCFT 
approximation. A linear constitutive relation was used: 
e((j>) = e + Aetp. Solid lines are zero-field phase bound- 
aries, while dotted lines marked with A, B, and C are the 
new phase boundaries in the presence of an external field. 
There appears a triple point where three phases meet. 



18 



At values of 4>q smaller than the triple point value, a di- 
rect transition between hexagonal and disordered phases 
occurs as a function of \ (or T). Part (b) is a phase 
diagram in the (x, Eq) plane for fixed 4>o- Dashed lines 
are a result of SCFT calculation, while the solid lines are 
from Eq. ( |37[ ) augmented with higher order terms in Ae 
and tfq. The higher-order corrections in Ae are essential 
since an interchange of the dielectric constant of the A 
and B polymers does not leave the electrostatic energy 
invariant. As the electric field is increased from zero, the 
area occupied by the distorted bcc phase shrinks in favor 
of the disordered and hexagonal phases. The phase lines 
terminate at a triple point (xtn E tI ). The bcc phase com- 
pletely disappears when E > E tr . The transition from 
spheres to cylinders was also studied for BCPs confined 
to thin films using SCFT. The deformation of the spheres 
from their perfect shape and the stability diagrams were 



calculated for the various system parameters (Matsen 
2006a|. 



In the above calculation, the zero-field critical point at 
(Nx,4>o) — (10.49, 1/2) is left unchanged because of the 
linear relation e{4>) and because lamellae and cylinders 
suffer no electrostatic penalty once oriented along the ex- 
ternal field. This is not true, however, if the constitutive 
relation has a nonvanishing quadratic dependence. To see 
that this dependence can change the ODT point, consider 
e" ^ in Eq. @, and ip = ip L cos(qx) [see Eq. Q], 



4>q = 1/2 (symmetric lamellae), and E is in the y or z di- 
rections. We find that / os = -(l/2)eJ5g - (l/8)e'V|£^. 
The electrostatic energy of the disordered phase is / cs = 
— (1/2)sEq, and therefore the electrostatic energy favors 
oriented lamellae over disorder if e" > 0; on the other 
hand, when e" < lamellae are expected to melt for 
large enough field even if they are oriented in the "right" 
direction. 

Lastly, we mention that the non-mean-field effect of 
uniform electric fields on composition fluctuations in 
symmetric block copolymers has been calculated as well. 
The weak first-order phase transition occurring in short 
polymers is shifted closer to the mean-field value, that is, 



the ODT point is shifted to higher temperatures ( 


Gunkel 


et al. 


2007 


Stepanow and Thurn-Albrecht 


2009 


1. As a 



result, phase separation is favored even if e = 0. The 
shift to the ODT temperature was estimated to be about 
2.5 K. A systematic measurement of the constitutive re- 
lation s((f>) is required in order to compare theory and 
experiment and to distinguish which of the terms - Ae 
leading to mixing, fluctuation effect leading to demixing, 
or e" leading to either mixing or demixing - is dominant 
in order-order phase transitions in polymer systems. 



V. LIQUID MIXTURES IN ELECTRIC FIELD 
GRADIENTS 



form fields, we defined uniform fields as those present 
in the system if the dielectric constant is uniform. In 
the language of Eq. (4]), this means E is constant. 
Spatially nonuniform fields have important consequences: 
they break the translational and the rotational symme- 
tries, and therefore the symmetry of the phase diagram 
around <p c is also broken. For a bilayer of two liquids in 
perpendicular electric field [Fig. [TJa)] a uniform field Eo 
deforms the interface only above a certain value, whereas 
a nonuniform field Eo deforms an interface no matter how 
small its strength is. Mathematically, when E is nonuni- 
form the linear terms in ip in the expansion of / cs [Eq. (|5|] 
survive the integration. The resulting diclcctrophoretic 
force tends to "suck" material with high e toward regions 
with high values of E 2 . 

As before, we stick for concreteness to a binary mix- 
ture of two simple liquids, denoted 1 and 2, with similar 
molecular volumes. The conclusions we draw are valid 
with only small modifications also to liquid-vapor coex- 
istence in field gradients. The mixture free energy as a 
function of composition, f m ((f>), is taken to be symmet- 
ric with respect to \<p — 1/2|. The mixture has an upper 
consolute point; if T > T c the energy is convex at all com- 
positions and the mixture is homogeneous. Below T c , f m 
has a double-minimum shape. The transition (binodal) 
curve (f> t {T) is given by df m (cj> t ,T) / d(j> = 0. 

We start with neat dielectric liquids, and in Sec. |V.B| 
we generalize to ion-containing liquids. 



A. Ion-free mixtures 

The mixture is subject to electric field emanating from 
a collection of conductors 1, 2, . . . , N, with prescribed 
potentials ipi, ■02, ■ ■ ■ > ipN, or charged with charges Q%, 
Qi, ■ ■ ■ , Qn each. We are interested in finding the 
potential distribution ip(r) and composition profile <ji>(r) 
as a function of the external potentials or charges. The 
total free energy density is given by 



/ = f m (4>) - \e{4>) (W>) 2 - ^ 



(38) 



The equilibrium profiles are governed by the following 
Euler-Lagrange equations: 

5f Sf m 1 de 2 n /o n x 

Si = ^-2M m) ~^ = °' (39) 



5£ 

Sip 



V (£(</>) Vi/O = 



In these equations, \x is the chemical potential needed to 
conserve the average mixture composition 4>q, 



(r)d 3 



(40) 



Up to this point we have been primarily interested in 
uniform fields. Since even inside an ideal parallel-plate 
condenser composition inhomogencities lead to nonuni- 



where V is the volume or is the reservoir chemical poten- 
tial in case of a grand-canonical ensemble. The second 
variational equation is simply the Laplace equation. The 
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set of differential equations [Eq. (39 1], are coupled and 
nonlinear and quite difficult to solve. Above T c , f m (4>) 
is convex. We therefore expect that under the spatially 
nonuniform forcing of the external field the composition 
</>(r) will be a smoothly varying function of space. Vari- 
ations in </> will be proportional to the variations in field, 
and if E (r) is continuous, so is <fi(r). 

Below T c , this is not true, and one can expect a differ- 
ent behavior as is shown below. Consider three "canon- 
ical" examples illustrated in Fig. [14] The first one is a 
charged isolated spherical colloid of radius R\ and total 
charge Q coupled to a reservoir at infinity. Azimuthally 
symmetric profiles depend on r only, where r is the dis- 
tance from the colloid's center. Thus cf) = <fi(r), e = e(r), 
and the field is given by E(r) = Q/ (4ire[<f>(r)]r 2 ) r . 
The second geometry is a charged wire of radius Ri 
and charge A per unit length, coupled to a reservoir at 
r — > oo. Alternatively, we may consider a closed con- 
denser made up of two concentric cylinders with fixed 
average mixture composition. In both cases the electric 
field is E(r) = A/ (2ne[4>(r)]r) f . A third system is the so- 
called "wedge," made up from two flat electrodes with an 
opening angle (3 between them, with potential difference 
V. In the wedge, E is given by E = V/((3r)6, irrespective 
of the profiles 4>(r) or e(r). The dielectrophoretic force 
existing in the wedge condenser may lead to discontinu- 
ous meniscus location in immiscible liquids (Jones 1974 1, 



or to exact cancellation of gravity effects in mixtures near 
their critical point ( Tsori and Lcibler 2005 1 . 



(a) 



OS) 





FIG. 14 Three model systems where field gradients lead to 
demixing. (a) A single charged colloid of charge Q and radius 
Ri . (b) A charged wire with charge density A per unit length 
and radius Ri or two concentric cylinders with radii -Ri and 
i?2- (c) A wedge comprised of two flat electrodes with an 
opening angle /3 and potential difference V. Ri and R2 are 
the minimal and maximal values of the distance r from the 
imaginary meeting point. 

Nonuniform electric fields will not demix liquids such 
that liquid interfaces are created parallel to the field. 
To see this, consider the two concentric cylinders with 
potential difference V. In the homogeneous state, e = 
e — 0o £ i + (1 ~~ < t ) o) e 2 is uniform. The demixed state 
has two uniform liquids of permittivities e\ and e 2 oc- 
cupying two portions of relative volumes </>o and 1 — (f>o, 
respectively, and the interfaces are parallel to E. The 
homogeneous and demixed states have exactly the same 
electric field and same value of F cs , but F m is higher for 
the demixed state. If the inner cylinder has a charge A 
per unit length immobilized to it, then the difference in 
.F es between the demixed and the homogeneous states is 
A 2 /(4 7 r)ln( J R 2 / J R 1 )[0 o /e 1 + (1 - o )/e 2 - 1/e], and this 



difference is always positive. 

The creation of interfaces perpendicular to the exter- 
nal field is in contrast to our intuition about uniform 
electric fields: despite the effects of dielectric interfaces, 
the dielectrophoretic force is the dominant mechanism. 
In addition, when a nonuniform field acts parallel to a 
liquid interface, it destabilizes it (Jones 1972 Perry and 



Jones 1975) 



Additional comment can be made on the possible mo- 
bility of charges inside the conductors: consider a single 
flat conductor at x = 0, charged with average charge a 
per unit area (uniform fields). The mixture is at the half 
space x > and has average dielectric constant e. Is 
it beneficial for the system to separate into two regions, 
one at y > 0, with surface charge a + 5a (at x = 0) 
and dielectric constant e + 5e (at x > 0), and one at 
y < 0, with surface charge a — 5a (y = 0) and dielectric 
constant e — 5s (x > 0)? Even without accounting for 
the reduced entropy of charges, this will not occur since 
the difference in F'es between demixed and homogeneous 
states is proportional to (5a/a — 5e/e) 2 and is always 
positive. Patches can be created on the conductor's sur- 
face, though, if the charges have effective repulsion due 



to van der Waals forces (Naydenov et al. 2007) 



Let us return to demixing in nonuniform fields. The 
nice thing about the three examples in Fig. [14] besides 
the geometrical simplicity, is the decoupling between the 
two equations [Eqs. (39)]. Essentially, we have already 
solved for Laplace's equation and only need to substitute 
the value of E in the variation with respect to (f>. We 
are left with a single nonlinear equation for <p(r), which 
has spatial dependence. Let us ignore for the moment 
the term proportional to (V</>) 2 in f m . This can be justi- 
fied because we are dealing with a volume term in E 2 (r). 
For a large system, the surface tension term acts only at 
sharp interfaces and can be neglected as a first approx- 
imation. For smaller systems, surface tension smooths 
out the composition profile <f>{r) and the neglect cannot 
be justified. 

In order to be concrete, consider the wedge geometry 
coupled to a reservoir at composition (j) (i? 2 ~~ * 00 ) an( i 
assume that e varies linearly with tfi W — in Eq.([3])]. 
We find the equation for ip is 



f' m {fp) = 



I* 



(41) 



(recall that ip — <j) — <p c is the deviation from critical 
composition). In the absence of field, the composition 
is (f>o , corresponding to a homogeneous phase above the 
binodal. With electric field, Eq. (41 1 gives an algebraic 
relation between ip and r. Above T c , the left-hand side 
is a monotonously increasing function behaving qualita- 
tively as p + ip 3, . Therefore, p increase monotonically as 
r decreases. This is the regular dielectrophoretic effect. 

However, below T c , the left-hand side has a sigmoidal 
shape; it behaves qualitatively as — <p + ip 3 . One may 
plot f' m {p) against ip. at a fixed temperature. The graph- 
ical solution to Eq. (41 ) is then given as an intersection 
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between this curve and a horizontal line independent of 
ip. At large values of r (small electric field), the intersec- 
tion occurs at negative values of ip. As r decreases, ip(r) 
increases. If the maximum field in the system V/(j3R\) 
is small, then the profile is still monotonic. If, however, 
V is large enough, the solution to ip to Eq. (41 1 occurs 



at negative values for large r's and positive values for 
small r's (large field). Necessarily, there is a location r 
where tp has a "jump." This is the field-induced phase 
transition: at a certain value of the voltage, the profile 
becomes discontinuous (Marcus and Tsori 2009 iTsori 



et at. 2004 1 . After demixing occurs, the wedge exhibits 



two regions with different compositions separated by a 
sharp interface. The two compositions of the coexisting 
phases are themselves nonuniform. 

The interface separating high and low values of (fi is 
at r = R (R\ < R < R 2 ). The critical potential V* is 
the minimal potential for demixing, and it occurs when 
R attains its smallest value, i.e., R — R\. As the po- 
tential increases, the interface moves to larger radii, and, 
in addition, the composition difference between the co- 
existing domains increases. For a closed wedge with 
average composition <fio, one readily finds that R(V) 
is bounded, and the maximum value of R is given by 
-^max = (1 ~ 0o)^i + 00-^2- If the wedge is coupled to 
an external reservoir and i?2 = oo, R grows indefinitely 
as a function of V. 

The dimensionless potential M w can be written as 
M w = V 2 NvqEq/ A[3kBT c R\. Here we generalize to mix- 
tures of polymers, N is the polymerization index and Nvq 
is the total chain volume. Close enough to the transition 
composition (fi t , the critical value of M w , M^, is bounded 
by the following expression: 



(fit - cf) T d 2 f m (<j> t ) 
4Ae T r d(fi 2 



g(R2/Ri) , 



(42) 



where we have used As = Ae/eq, and f m — Nvof m /kBT. 
g(x) = 2(x 2 — l)/(x 2 — 1 — 21ns) is a geometry-specific 
dimensionless function of the ratio x = R^/Rx- Fig- 
ure 15 ^a) shows M£j as a function of difference between 
the transition composition <fi t and 4>q. Solid lines are 
Eq. (42 1, while symbols are taken from a more accurate 
numerical calculation. Clearly, M w increases almost lin- 
early with <fi t — (pQ. The slope decreases at T approaches 
T c . Note that the similarity between the wedge capaci- 
tor and a rapidly rotating centrifuge allows us to predict 
phase-separation transitions for mixtures in centrifuges 
if the rotation frequency is higher than a certain critical 
frequency ( |Tsori and Leibler] |2007b| | . 

The behavior in the cylindrical (2D) and spherical 
(3D) cases is similar to the one described above for 
the wedge. There are expressions for the critical val- 
ues of the dimensionless charges for the cylindrical and 
spherical symmetries: M c = A 2 Nvo/16irkBT c R 2 £o and 
<5 2 -/Vwo/647rfc B r c i?fe 0l respectively. Figure 15 'b) 

as 



M, 

shows the critical values for demixing, M* and M, 
a function of temperature for three different composi- 
tions. Open symbols are calculated values of M* for 



two concentric cylinders, whereas filled symbols are M* 
for an isolated spherical colloid. For a given composi- 
tion, curves only start above the transition temperature 
(T/T c ~ 0.865, 0.944, and 0.986 for (fi = 0.2, 0.3, and 
0.4, respectively). Note that in general, M* and M* are 
much larger than - this is again because of the di- 
electric interfaces: in the spherical and cylindrical sym- 
metries, E is parallel to V</>, and therefore the dielec- 
trophoretic force oc As has to be large enough in order 
to overcome the penalty associated with dielectric inter- 
faces, oc (Ae) 2 . In the wedge E _L V<f>, and is corre- 
spondingly smaller. 
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FIG. 15 Critical demixing fields, (a) Scaled critical voltage 
for demixing for a closed wedge as a function of distance 
from the transition (binodal) composition (fit. Symbols are 
numerical values and line is approximate analytical expres- 
sion [Eq. j42|]. (b) Scaled critical charge of spherical and 
cylindrical colloids as a function of temperature. Open sym- 
bols are M* for a closed cylindrical system of compositions 
4>o — 0.2 (squares), (fio = 0.3 (circles), and (fio = 0.4 (crosses). 
Filled symbols are for an isolated spherical colloid cou- 
pled to a reservoir with the same compositions. Reprinted 

( |2008| ). Copyright 2008, 



with permission from Marcus et al. 



American Institute of Physics. 

How much should a 1 /im colloid be charged to cre- 
ate a demixing layer around it? We take the molecular 
volume wo to be 3 x 10~ 27 m 3 , N — 1 (simple liquids), 
and T c = 350 K. We find that M* = 0.001 leads to 
Q ~ 10 4 e. This may seem like a large number; how- 
ever, if the mixture contains even a low molecular weight 
polymer, then N = 500 gives Q ~ 470e. The enrichment 
layer created around the colloid may be quite important 
for the thermodynamic behavior of suspensions because 
field-induced demixing leads to capillary attraction, pro- 
moting colloidal aggregation. 

Finally, for a given value of (fio , there is an approximate 
expression for the range of temperatures AT above the 
transition temperature T t ((fi ) for which a given voltage 
is able to demix, 



AT: 



Nv 
2k B 



Ae 



V 



(43) 



This formula, applicable to the wedge, is valid close 
enough to the binodal. It should be contrasted with the 
Landau expression [Eq. ([36])] . When the factor N/2 is 
put aside the two formulas have similar form: e"Eq in Eq. 
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(36 1 is replaced by Ae(V//3Ri) 2 , and both expressions 
have similar magnitudes since Ae/e" ~ 1. However, the 
denominator of Eq. ( 43 1 has the difference between </> 



and the critical composition. This factor means that AT 
here is at least twice as large as in the Landau case and 
typically is much larger: \<p c — </>o| _1 ~ 20 or even more. 

The phenomenon described here has many similarities 
to regular prewetting. Important differences are (i) E is 
determined from Laplace's equation and a similar equa- 
tion does not exist for the case of long-range van der 
Waals forces between liquids and solid surfaces, (ii) In 
addition, the electric field needs not attain its highest 
value at the surface. Thus, "electroprewetting" is a spe- 
cial case of field-induced demixing. Consider an analo- 
gous system of a body of nonuniform heat conductivity 
k(t) enclosed by a surface £ with prescribed tempera- 
ture Tx;. Temperature is the analog of ip and k is the 
analogue of e. A question follows: is it possible to have a 
distribution k(t) such that the steady-state temperature 
profile T(r) achieves a maximum or a minimum at the 
bulk rather than at E? A positive answer would mean 
that a confined mixture could spontaneously demix due 
to heat flow even without Soret or other effects. The 
answer to the question is negative. Since T is governed 
by V(kVT) = 0, we find Vk • VT + nV 2 T = 0, and ex- 
trema do not satisfy this equation because VT = and 
V 2 T 7^ 0. However, the point here is that electric fields 
are different: the fact that if) is largest on the electrodes 
does not prevent E 2 from being the largest elsewhere. 

Up to this point, the liquids were considered as pure 
dielectrics. In the following section the phase-separation 
phenomenon is described in the presence of ions. 



B. Phase separation in ion-containing liquids 

The field-induced phase transition described above 
originates from electric-field gradients. In nonpolar 
liquids, such gradients occur when the electrodes are 
curved. In one dimension, e.g., in parallel-plate con- 
denser, field gradients are forbidden by the equation 
V • E = 0. For curved electrodes, the geometry dictates 
the length scale over which the electric field falls off. On 
the other hand, mixtures of polar liquids (e.g., aqueous 
solutions) contain some amount of charge carriers. In 
such mixtures, the physics is rich and quite different from 
the simple dielectric case. The most important feature is 
due to screening, brought about by dissociated ions ac- 
cumulating at the electrodes. Screening means that the 
electric field is substantial only close to the electrodes, 
within the screening distance A. Field gradients origi- 
nate from both geometry and screening, and the phase 
transition depends on two lengths. The ionic screening 
thus adds to the dielectrophoretic force which separates 
the liquid components from each other. Since screening 
is omnipresent, phase separation may occur even near 
parallel and flat charged surfaces, i.e., in one dimension. 

But ions have another effect besides increasing the di- 



electrophoretic force. Ions have in general different sol- 
ubilities in the different liquids. As an ion drifts toward 
the electrode, it might "drag" with it the pre ferred liq- 
uid component (Onuki and Kitamura 20041. Usually 



the two ionic species are preferential in the same liquid, 
but it is also possible that they are preferentially soluble 
in different liquids. The solubility introduces a force of 
electrophoretic origin, proportional to the ions' charge. 
In binary mixtures, the parameter Au + measures how 
much a positive ion prefers to be in the environment of 
liquid 2 over that of liquid 1. To lowest order, the inter- 
action per unit volume of positive ions and the mixture 
is given by — Au + n + (r)(f>(r) , where n + is the positive ion 
number density. A similar expression exists for the neg- 
ative ions. The system free energy density including the 
interaction of ions and field is given on the mean-field 



level by (Onuki 2006 Onuki and Okamoto 2009) 



/ = 



fmtt) " \<<t>) (VV>) 2 + (n+ - n-) eif> (44) 



+ ksT \n + In (vqii + ) + n In (v^n )] — \ + n + 



A n 



[i(f> — (Au + n + + Am n 



const. 



In the above, A and /i are the Lagrange multipliers 
(chemical potentials) of the positive and negative ions 
and liquid concentration, respectively. 

The free energy needs to be extremized with respect to 
the fields cf>, ip, and in keeping with a fixed ion concen- 
tration V -1 J n ± (r)d 3 r = no. The Poisson-Boltzmann 
equation is obtained from the variation in / with respect 
to the potential and is coupled to Sf /8(f). This formalism 
is similar to other formulations, but the crucial difference 
is the bistability of the mixture energy f m (Ben-Yaakov 



et al. 2009 Onuki and Kitamura 20041. 



Solution to the problem is difficult to obtain because 
of the essential nonlinearity: a linear Poisson-Boltzmann 
equation does not suffice since imposing eif> -C fceT 
means small potential and thus forbids the large elec- 
tric fields necessary to drive the phase-separation tran- 
sition. The transition temperature T t (cf>o), defined by 
dfm(<t>o,Tt) / d(f> — 0, is shifted by an amount AT due 
to the external potential. An approximate formula for 
AT near a single planar electrode at x = can be ob- 
tained by the assumption that field gradients are domi- 
nated by ionic screening rather than by e variations. It 
reads (Tsori and Leibler 2007a I 



AT 



lAd 



Au 

knT r 



Nn v 



\4> c — (f>o\ P \ k B T r 



eV 



, (45) 



where Au = Au + — —Au~ and V is the electrode poten- 
tial. This expression holds as long as T t + AT is smaller 
than T c ; at all temperatures T > T c , the composition pro- 
file <f>{x) varies smoothly with no abrupt jump. The shift 
in T t should be compared to Eq. ( 43 ) for a mixture with- 
out ions and to Eq. (36 1 for a mixture in uniform field. 



The dielectrophoretic and electrophoretic driving forces 
appear to have similar significance since both Ae/s c and 
Au/UbT are ~ O(l). The numerator has a small factor 
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NtioVq; in pure water and if we take Nvo = 3 x 10~ 27 
m 3 , we get Nnovo ~ 10~ 7 . Small as this factor is, and 
even ignoring the denominator of \<p c — 4>q\, the exponen- 
tial factor exp(eVyfcsT) is the dominant term in AT. It 
is usually huge, recalling that at room temperature 1 eV 
equals 40k B T. 

One can therefore expect that in almost all circum- 
stances the electroprewetting discussed above should oc- 
cur around a charged colloid in aqueous solutions. The 
width of wetting layer is small but may nonetheless influ- 
ence the effective interaction potential between two col- 
loids orbetween a colloid and a nearby surface (Bauer 



et al. 2000 1 and the thermodynamic behavior of suspen- 



sions of charged particles (Beysens and Esteve 19851 



Electric field gradients may also play a role in such phe- 
nomena as the "critical Casimir effect," recently reported 



experimentally by Hertlein et al. ( 2008 1 



The above wetting transition from a homogeneous to a 
demixed solution has several aspects. First, if the two or 
more constituent pure components have sufficiently dif- 
fering indices of refraction, the optical properties of the 
solution will be drastically different in the two states. 
The demixing of an initially transparent solution leads to 
the creation of optical interfaces, and these can be used 
to scatter, reflect, and refract light. Particularly in a 
microfluidic environment, the phase separation could be 
utiliz ed to create wavegui des which are liquid and switch- 
able (Psaltis et al. 20061. The mechanical properties of 



the mixture are also different whether it is homogeneous 
or not. If the two pure components have different vis- 



cosities, as is the case of the experiments (Tsori et al. 



2004), then after demixing a thin lubricating layer is cre- 



ated at the charged surfaces, and this thin layer com- 
pletely alters the flow profile between the surfaces. Fi- 
nally, the field-induced phase separation can be used to 
control chemical reactions of, say, two molecules A and 
B reacting in a mixture of liquids 1 and 2. The depen- 
dence of concentration on electric field can accelerate the 
reaction or conversely stop it. The spatial dependence 
can be used to confine the reaction to small volumes or 



to two-dimensional interfaces (Tsori and Lcibler 2007a). 



C. Phase separation in surfactant mixtures 

Phase transitions in nonuniform electric fields are not 
restricted to demixing in liquid mixtures or to a liquid- 
vapor coexistence. In a set of beautiful experiments, Lee 
et al. ( |1994 | demonstrated that a flat monolayer of a 
binary lipid mixture on the water-air interface undergoes 
a transition from a homogeneous to a phase-separated 
state when an insulated wire passing perpendicular to the 
monolayer is charged. The spatially varying electric field 
couples to the different dipole moments of the constituent 
molecules and attracts the molecules with large dipole 
toward the wire. 

These experiments have an important difference from 
the demixing described before: the electric field couples 



to the fixed dipole of the lipids. Thus, the effect scales 
linearly with E rather than quadratically. Consequently, 
when the wire potential is V, lipid 1 is attracted to the 
wire, whereas the potential reversal V — > —V reverses the 
effect and accumulation of lipid 2 is observed at the wire. 
Note that the dipole orientation is restricted to be in one 
direction; if dipoles could flip upside down, a reversal 
of the field E — > —E would leave the system invariant 
and the effect would again scale as E 2 . Similarly to the 
demixing in simple liquids, below the critical pressure 
(which replaces T in the phase diagram) only smooth 
gradients in <p appear, whereas below the critical pressure 
true demixing occurs, as evidenced by the sharp interface 
separating coexisting domains. 



VI. CONCLUSIONS 

In this Colloquium we considered liquids and polymers 
in external electric fields. The key quantities for the phys- 
ical behavior are the linear and quadratic coefficients in 
the expansion of e(<j)), namely, Ae, and e" in Eq. 
and the conductivity a. We discussed the significance 
of these quantities in several phenomena occurring when 
the field is either spatially uniform or nonuniform. We 
have emphasized that in a uniform field Eo, there is a 
preference of electrostatic origin to have dielectric inter- 
faces parallel to E . This preference, depending on Ae, 
gives rise to the perpendicular field instability observed 
for liquid films discussed in Sec. [n] The wavelength of 
instability A* is the fastest-growing mode of a dynam- 
ical evolution. What is the smallest feature size possi- 
ble in these experiments? The answer is not clear, and 
an effort is undertaken by several groups to minimize 
A* by optimizing the fluid's surface tension, film thick- 
ness, and other parameters. Dielectric interfaces are also 
responsible for block-copolymer orientation and order- 
order phase transitions in electric fields. Current effort 
focuses on finding ways to eliminate in-plane defects and 
enlarging the grain size. Destabilization of copolymer in- 
terfaces and polymer brushes may also occur, but here 
the wavelength is dictated by equilibrium considerations. 
In addition, a uniform electric field deforms the shape of 
suspended drops and of drops on solid substrates. 

We have outlined some of the effects related to disso- 
ciated ions, commonly present in most simple and poly- 
meric liquids. Finite conductivity introduces additional 
stress and can lead to oblate shapes of an isolated drop. 
Electrowetting is a result of charge accumulation at the 
droplet-substrate interface. The reduction in the droplet- 
substrate interfacial tension leads to an apparent change 
in the contact angle, and this has been widely used in 
applications. Ions also play an important role in copoly- 
mer orientation and phase transitions. Complexed ions 
mainly contribute to enhancement of Ae, whereas mo- 
bile ones may give rise to wettinglike effects due to their 
preferential solubilities. An experimental effort is under- 
way, aiming to understand when complexed ions are more 
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dominant than mobile ones and vice versa. The study of 
electric-field-induced wetting of complex phases, such as 
block copolymers, is still at early stages and deserves ex- 
perimental and theoretical treatments. Ions introduce 
stresses when the fields are time varying, and the result- 
ing torques tend to orient the copolymer domains. The 
dynamics of these out-of-equilibrium systems is an open 
problem that also requires further investigation. 

Uniform electric fields can also lead to phase transi- 
tions in liquid mixtures by the Landau mechanism. These 
transitions are insensitive to Ae but rely on nonvanishing 
e" . Depending on the sign of e", the whole binodal curve 
can be shifted to higher or lower temperatures, but this 
shift is rather small, ~ 0.03 K. 

Spatially nonuniform fields deform ordered block- 
copolymer phases due to a nonvanishing Ae. The copoly- 
mer structure, a compromise between elastic bending and 
compression energies and the electrostatic energy, has 
not received enough attention so far. It is hoped that 
research in this direction can show us how to eliminate 
defects and perhaps even how to tailor lamellae or other 
phases in a prescribed way on a substrate. In the much 
simpler systems of liquid mixtures, bending is irrelevant, 
and nonuniform fields have a strong effect on the phase 
behavior. The direct coupling between field and composi- 
tion leads to a phase transition which depends on Ae and 
cannot be simply described by a renormalization of T c . 
The change in transition temperature is much larger than 
in the Landau mechanism. Wetting layers are created 
around curved charged objects but can also appear far 
from a solid interface. The presence of dissociated ions 
greatly enhances the transition and allows it to happen 
at flat objects. This transition can have numerous appli- 
cations in microfluidics, in micro and nanorhcology and 
lubrication, in chemical reactions, etc. The current study 
is only at its infancy and there is still plenty to learn: for 
example, what is the size and velocity of liquid domains 
appearing upon switching the field on, why interfacial in- 
stabilities are observed after demixing takes place, what 
is the role of critical fluctuations and the non-mean-field 
behavior, how does heating affects the transition, etc. 

The phenomena outlined above are fundamentally in- 
teresting and also have numerous applications. We are 
certainly going to learn more about this active field of 
research in the years to come. 
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